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Abstract 



Lattice paths effectively model phenomena in chemistry, physics and probability theory. 
Asymptotic enumeration of lattice paths is linked with entropy in the physical systems 
being modeled. Lattice paths restricted to different regions of the plane are well suited to 
a functional equation approach for exact and asymptotic enumeration. This thesis surveys 
results on lattice paths under various restrictions, with an emphasis on lattice paths in the 
quarter plane. For these paths, we develop an original systematic combinatorial approach 
providing direct access to the exponential growth factors of the asymptotic expressions. 
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Background 



Chapter 1 



Introduction 



In many applications, lattice paths are models of physical and chemical phenomena, with 
the enumerative results closely linked with physical properties. For statistical mechanical 
systems modeled using combinatorial structures, growth constants, known as exponential 
growth factors in the enumerative world, provide information about the entropy in the phys- 
ical system IfTSll . 

The focus of this thesis is a family of planar lattice walks confined to different regions 
of the plane, with an emphasis on the quarter plane, a region for which explicit closed form 
enumeration is presently unknown for many members of this family. For many models of 
quarter plane walks, asymptotic expressions are known, but the present methods of finding 
them require heavy computation. 

Motivated by finding a more intuitive strategy for proving the exponential growth factors 
combinatorially, we study planar lattice paths with small or nearest neighbour steps. These 
are paths that move from the origin to another point in a region R c Z^, using a set of 
steps S taken as a subset of S = {-1, 0, 1}^ \ {(0, 0)}, which we also may represent by the 
corresponding compass directions 

S = {N, NE, E, SE, S, SW, W, NW}. 

The number of steps taken by a walk is the length of the walk, and we count the number 
of walks of a given length on a step set subject to some restrictions. Figure II. H is an 
unrestricted lattice path of length 1000 with steps taken from {N, E, S, W}, starting at the 
origin. 
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results. 



CHAPTER 1. INTRODUCTION 
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Figure 1.1: A walk of length 1000 on S = {N, E, S, W}. 



Our approach uses generating functions, that is, formal power series of the form 

Rs(t) = J]rs(n)t", 

n>0 

where r3(n) is the number of walks of length n taken with steps from S in the desired 
region. We index the generating functions by the step set the lattice paths are taken on, 
and the region that the paths are restricted to will be clear from the context. The growth 
constants, or exponential growth factors, are found in the asymptotic expression for the 
sequence (re(«)). For the cases considered here, these expressions take the form 

rs(n) ~ K/3"e(n), 

Where a: is a positive constant, /3 is the exponential growth factor and 6(n) is some subex- 
ponential growth factor. We define this formally later, but in an intuitive the asymptotic 
expression is an approximation of the sequence re(«), the quality of which increases with 
n. The form taken by 9{n) can vary, and is related to the kinds of equations that i?s(0 
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satisfies. We focus on generating functions classified as rational, algebraic, transcendental 
D-finite and non-D-finite. For some classes of generating functions, the form that 6{n) takes 
is known: for example, for algebraic functions 6{n) is a Laurent monomial, and for rational 
functions k = 6{n) = 1 . We give a summary of the current state of research, including our 
contribution, of the four classes of planar lattice paths discussed in this thesis in Table ll.ll 

In fj\, Bousquet-Melou and Mishna show that there are 79 non-equivalent models of 
walks confined to the quarter plane. They then classified these according to the cardinal- 
ity of a group associated with the step set. For 23 models, the groups are finite and the 
associated generating functions are proven to be D-finite, that is, they satisfy a linear differ- 
ential equation. Unfortunately the proof does not give the differential equation explicitly, 
and even if it did, there would be no guarantee of finding a solution. They conjectured the 
remaining 56 to be non D-finite, and this was very recently proved for 51 cases by Kurkova 
and Raschel in [|H. 

We provide here a combinatorial method which will confirm the exponential growth 
factors numerically obtained by Bostan and Kauers in p1| for all of the 23 D-finite cases. 
The data provided by Bostan and Kauers is an invaluable guide (as it usually is when you 
know what you're trying to prove), but often the exponential growth factors can also be 
hypothesised using simpler numerical methods. They show in [4] that these generating 
functions are not only D-finite, but are also G-series. This means that the number of walks 
grows like 

K/S"}!" (log ny 

for some K,a,/3,y. By generating the sequence using experimental techniques and taking 
the ratio of successive terms far along the sequence, we can approximate /3. For instance, 
if we are looking at a counting sequence with an algebraic generating function, the ratio of 
two successive terms would be 

win + l) K/3"^\n+\y _ Jn + lV 
w(n) K/3"n" \ n ) 

Using this, we can get an approximation foryS that decreases in quality as a grows. In [4], 
Bostan and Kauers never encounter an exponent a larger than -1/2, so taking a = I would 
provide a reasonable worst case scenario. For n = 1000 terms, this gives 1.001/3 as a result, 
and for n = 2000, 1 .0005/3. In general, this is not precise enough that a computer could 
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guess the growth constants, but if the result is close to the growth constant of a related 
model, then it seems clear to a human what these are converging on. 

The method for the combinatorial proof of these 23 factors appeals to relaxations and 
restrictions of the models. By considering the right relaxation (respectively restriction), 
one can place an upper (resp. lower) bound on the exponential growth factor. Elementary 
analysis does the rest. 

In Chapter |2l we introduce the formal machinery of lattice paths, and some supporting 
notation that we find useful throughout the text. These are at the most general level only, 
and more specific definitions will be given in later chapters. 

Chapter |3] will give an outline of the link between exponential growth constants and 
statistical mechanical properties. We also discuss the famous Self Avoiding Walk (SAW) 
model in three different configurations, and the current state of knowledge on bounds for 
the empirically known exponential growth factor. 

In ChapterlH the necessary elements of analytic combinatorics are surveyed. This starts 
from an introduction of a symbolic method, and ends with a proof of the correspondence 
between the radius of convergence of a generating function and the exponential growth 
factor of the coefficient sequence. 

Chapter |5] reviews planar lattice walks with small steps which are completely unre- 
stricted. This is the classical case, easily enumerated with a rational generating function. 

Chapter |6] is the first restriction to the general class, where we consider directed paths. 
These paths have a priveleged direction of increase, with the x-coordinate of the endpoint 
equal to the length of the path. We survey results from [l2], with an emphasis on directed 
paths confined to the upper half plane. Generating functions are shown to be algebraic 
through the use of the kernel method, and asymptotic enumeration is given. This chapter is 
significant, as it contains the kind of result that we are working towards: exact enumerative 
formulas dependent on the drift parameter. 

Chapter |7] considers the class of walks confined to the upper half plane. A bijection 
between this class and the class of directed paths from Chapter |6] is given. This shows 
generating functions to be algebraic, and enumerative results follow immediately from the 
results of flU given in Chapter |6] 

Chapter [8] covers the main problem, planar lattice walks with small steps confined to 
the quarter plane. The class is split into two families, those known to be D-finite and those 
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known or conjectured to be non D-finite. Here we give an original, systematic method of 
proving the exponential growth factors for 23 D-finite models of quarter plane lattice paths. 
Chapter |9] concludes the thesis. We consider the generalisations of the technique to all 
non-equivalent quarter plane models. We also consider the prediction of the polynomial 
growth of the counting sequence, and how we can put this work into a single, general 
framework. 



Chapter 2 



Definitions and Notation 



This chapter introduces formal definitions relating to planar lattice paths with small steps. 

Definition 1. Fix a finite set of vectors in Z^, S = {{ai,bi), {a2, ^2), •••, {cim, bm)}. A lattice 
path or walk taken on S is a sequence v = (vi, V2, ..., v„) such that each v,- is in S. We may 
realise this geometrically as the sequence of points {po, pi, ..., p„) such that po = (0, 0) and 
the vector difference p, - p,_i = Vj. The elements of S are known as the steps, and the 
number of steps in the sequence v, n, is called the length of the path. 

We recall the definiton of a Laurent polynomial which encodes the step set S in a useful 
way. 

Definition 2. For a set of steps S = {{ai,bi), {a2, ^2), •••, {cim, bm)} we define the inventory, 
or characteristic polynomial, of S to be the Laurent polynomial denoted by S {x, y) which 
is defined by 

m 

S{x,y):=Y,^9'- 

(=1 

Sometimes it is useful to add weights to the steps of a walk. In combinatorial applica- 
tions, these weights are often natural numbers. Here they correspond to different coloured 
steps in the same direction. In probabilistic applications, the the total weight of a step set 
(the sum of weights over the steps) would be equal to 1 . In either case, we define a set of 
weights n = {w\, ...,Wm}, and modify the inventory to include the weights: 



^(•^,j) = £>v.-^V'. 



1=1 
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In order to simplify analysis we restrict ourselves to what are known as small steps 
until the end of Chapter [8] This is just a way of saying that at each step, the walk moves 
only to nearest integral lattice points diagonally, horizontally and vertically. We discuss the 
generalisation to larger steps in Chapter |9] 

Definition 3. If S is a subset of S = {-1, 0, 1}^ \ {(0, 0)} then we say that S is a set of small 
steps. 

Sometimes, rather than using a vector representation it is convenient to represent a step 
set pictorially, as in Figure IZTl or by the corresponding points on the compass. So, moving 
clockwise from the (0, 1) step in Figure O these would be N, NE, E, SE, S, SW, W, NW, 
with corresponding inventory 

X I \ I y 

S {x,y) = y + xy + X + - + - + — + - + -. 

y y xy X X 



^ 



Figure 2.1: The full set of small steps S. 

We can also represent a walk with the corresponding polygonal line in Z^. Shown in 
Figure 12.21 is a geometric representation of a typical unrestricted walk on the set of small 
steps S = {NE, SE, SW, NW}. By unrestricted we mean that we allow self intersections 
and the final point pn to be anywhere. 

We find it useful at times to take the partial sum of a Laurent series or polynomial in 
which the exponents of a variable are negative. We recall some notation for this. 

Definition 4. Let H(y) be a Laurent series or polynomial in y with coefficients taken 
from Z{x). We denote by {y^^}H(y) the projection onto the pole part of H{y): the partial 
sum of H{y) where all terms contain a negative index of y. 

We illustrate this by example after the following definition. 

Definition 5. We say that a step {i, j) is ;c-positive if / > 0. Similarly, (/, j) is ;c-negative 
if i < 0, y-positive if j > and y-negative if j < 0. 
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S = X, « = 7 



Figure 2.2: A walk of length 7 on a typical set of small steps S. 



Example 2.1. Take S (x, y) to be the inventory of the complete set of small steps. Then 



{y''"}S(x,y) = {/"} \y + xy + x + - + - + — + - + - 

y y xy X X 



X 1 1 

- + - + — 

y y xy 



is the Laurent polynomial representing all of the y-negative steps. 



Finally, we formally define what we mean by asymptotic growth. As mentioned in the 
introduction, intuitively it is a continuous function of n which is an approximation of a 
sequence (a„), the quality of which improves as n increases. 

Definition 6. Let a : N ^ N be a function (such as the coefficient sequence of a generating 
function associated with a combinatorial class). We say that a„ grows like or has asymptotic 
growth g{n), g -.W^ ^ W^, denoted 



a(n) ~ g(n) 



if 



lim 



= 1. 



This concludes the general definitions for planar lattice paths with small steps. In Chap- 
terlH more definitions are given as the formal framework of analytic combinatorics is built, 
and in Part [III more specific definitions will be made for each class, as restrictions are 
introduced and more machinery becomes necessary. 



Chapter 3 



Statistical Mechanics 



Recent activity in lattice path combinatorics (for example, [I3l|4l[5l|71|22l|) is due in part to 
the success with which they are applied to a wide variety of areas. Two often cited examples 
are applications to other combinatorial models through bijections, such as between lattice 
paths and plane partitions in the work of Alegri et al. in HJ, and queueing theory, as in 
Bohm's work [Q on lattice path counting and the theory of queues, but the example which 
will form the main focus of this chapter is the application of lattice models to the modelling 
of physical phenomena through statistical mechanics. 

These applications take different forms, and a nice introduction to a number of these 
can be found in Janse van Rensburg's book 'The statistical mechanics of interacting walks, 
polygons, animals and vesicles'" filSil . Of these applications, the most significant to us is the 
self avoiding walk (SAW). These walks are defined in the same way as our walks in Chapter 
[21 but we make the restriction that in the sequence of points {po, ..., p„) there is no i + j such 
that Pi = pj. That is, the walk may not visit the same place twice. This makes these walks 
a good model for linear polymers, something we explain more in the next section. First, 
we should explain why combinatorial models are well suited to physical problems such as 
these. 

Information about the entropy of a physical system is available in the form of a model's 
enumerative information. In statistical mechanics, entropy is a measure of the number of 
ways in which a system may be arranged, often taken to be a measure of 'disorder' (the 
higher the entropy, the higher the disorder). This definition describes entropy as being 
proportional to the natural logarithm of the number of possible configurations of the indi- 

11 
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vidual atoms (with Boltzmann's constant as the proportionahty constant). In order to find 
the number of possible configurations, we would be forced to count them. This is where 
enumerative combinatorics comes in. Since these are purely mathematical models we ig- 
nore the proportionality constant and any units of measurement, and focus solely on the 
enumeration. 

So, if S is a step set, and f is a class of lattice paths on S modeling some physical 
phenomena, we can find the combinatorial version of the entropy E(n) by taking 

E(n) = log p(n), 

where p{n) is the number of lattice paths in f of length n. However, sometimes in statistical 
mechanics it is desirable to understand the behaviour of a model as the size approaches 
infinity. The sequence p{n) is unbounded in general, so if this limit were taken the entropy 
would also be unbounded. There is another measurement though, related to entropy, called 
the limiting free energy or the thermodynamic limit. It is given as 

F = lim -E{n). 

n—>oo fl 

Now for many lattice paths (including those found in this thesis) the counting sequence has 
growth 

p{n) ~ ra'^yS", 

for K,/3 eW^ and or < 0, so the limiting free energy would be 

F = lim - logim^fi") = /3, 

n—>oo fl 

which is also known as the exponential growth factor of the number of states. So by finding 
exponential growth factors for models, we can give information about the entropy of a 
statistical mechanical system. 

In the following section we survey a few different self-avoiding lattice paths and their 
applications. In light of the link between entropy and exponential growth, we also give the 
asymptotic expressions for each model. 

Self-avoiding walks 

A polymer is a connected arrangement (similar to a graph-theoretic structure) of molecules 
known as monomers. Polymers can take on tree-like structures, have internal polygonal 
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Figure 3.1: A SAW which cannot be extended to form a longer SAW. 



structures, or simply be a linear chain of monomers. A linear polymer is a polymer where 2 
monomers have only one neighbour and the remaining monomers have exactly 2 neigh- 
bours. As already mentioned, self avoiding walks (SAW) make very succesful models of 
linear polymers suspended in solution, see [|23l[T6l l9l for some examples of this application 
in 2 and 3 dimensions. We focus on 2 dimensions here. 

We call the class of self avoiding walks Ji, and they are defined as in the beginning of 
this chapter. There is a problem though: this class isn't amenable to recursion. That is, 
not all SAWs can be added to in order to create a SAW of longer length. Figure ISTI shows 
an example of a SAW on S = {N, E, S, W} which cannot be extended to create a SAW of 
longer length: any extension of the path will violate the self avoidance restriction. A lack 
of a recursive structure means that we cannot build a specification and hence a functional 
equation satisfied by a generating function, as is done in Chapter |4] Thus, many of the 
enumerative results on SAW are empirical, including the asymptotic expression for the 
number of walks of length n: 

a(n) ~ /^(2.638f n"^^^ 

where /c is a positive constant [|T4l . 

If we make restrictions on the class, simplified models can be created which either are 
automatically self avoiding (directed paths) or are able to be attacked recursively with self 
avoidance enforced (partially directed paths). These simplifications allow researchers to 
work on self avoiding combinatorial models in the usual way, finding rigourous results on 
related problems and forging new techniques for attacking an old, stubborn problem. 

Even very simple combinatorial models have been succesfuUy applied in this way. One 
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example is the application of directed paths to the statistical mechanics of linear polymers 
subject to a force, such as in an atomic force microscopy (AFM) experiment. We show an 
example of a directed path in Figure 13.21 Atomic force microscopes are able to interfere 
with a substance at an atomic scale, so in linear polymer experiments they would dip into 
a solution and hold on to a polymer by a single monomer, applying a force in a vertical 
or horizontal direction, making the polymer take a shape similar to one of these types 
of paths. The methodology is to analyse combinatorial models for phase transition and 
entropic information through asymptotic analysis of associated generating functions, and 
confirm with AFM experiments 



> 



Figure 3.2: A directed path on S = {NE, E, SE}. 



We also find rigorous bounds on the exponential growth of SAW by considering a sub- 
class. Take the class P of partially directed paths on the square lattice. This is a walk taken 
on the set S = {N, E, S} with self avoidance enforced. An example is shown in Figure |331 

The walks in P form a subclass of general SAW on the square lattice, and can be solved 
explicitly, with generating function 

t(l-t) 



Pit) 



\-2t-f 



The dominant singularity for this generating function is a simple pole at V2 - 1 . As will be 
explained in ChapterlH we can use singularity analysis to compute 



p{n) ~ K 



V2-I 



k(2AUT. 
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Figure 3.3: A partially directed path provides a lower bound on the exponential growth of 
SAW. 



Using the techniques we discuss throughout this thesis, we can use this asymptotic estimate 
to provide a rigorous lower bound on the exponential growth factor for all SAW on the 
square lattice. 

Combinatorics could be used to analyse other thermodynamic properties of a system. 
At any moment, a system is distributed across an ensemble of A^ microstates, each denoted 
by / and having a probability of occupation p,, and an energy Ei. These microstates form 
a discrete set as defined by quantum statistical mechanics, and Ei is an energy level of the 
system. One may then look at the intemal energy of a system, which is the mean of the 
energies Ei, or the heat and work of a system. Moreover, lattice statistical mechanics is 
not restricted to systems modeled by lattice paths [[TSl . This makes the interface between 
combinatorics and statistical mechanics a rich junction, with combinatorial approaches not 
limited to enumeration. 



Chapter 4 

Analytic combinatorics 



When trying to enumerate a family of objects which is amenable to recursion, like lattice 
paths, a standard tool is the generating function. These were originally introduced as formal 
power series, or as Wilf colourfuUy puts it in ^2E\, a 'clothesline to hang a sequence upon,' 
these tools have proved very useful in many applications. Stanley has written two volumes 
Il24l |25]| on the topic of enumerative combinatorics in general, covering applications of 
generating functions to algebraic structures as well as the algebraic nature of the objects 
themselves, providing in [|24]| a classification of formal power series which we shall employ 
a little later on. 

In IfTTTl . Flajolet and Sedgewick take the field much further, considering the formal 
power series inside their radii of convergence as complex analytic functions: no longer 
just algebraically. This allows the complex analysis theorems of Cauchy to be applied 
with impressive enumerative results, obtaining asymptotic expressions for the coefficient 
sequences of the generating function with little effort (once the theorems have been proved). 
They describe their approach as an operational calculus for combinatorics organised around 
three components: symbolic methods, complex asymptotics and random structures. The 
following gives an introduction to the symbolic methods and complex asymptotics, drawing 
on elements from [|Tni28ll25l. 



16 
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4.1 Combinatorial classes and ordinary generating func- 
tions 

We follow the exposition of the symbolic methods given in ifTTI . concentrating on the points 
useful to us. The first of these will be the combinatorial class. Earlier, we mentioned that we 
aim to enumerate, either exactly or asymptotically, subclasses of 'Wq, the class of planar 
lattice walks on a given step set S. 

Definition 7. A combinatorial class or simply class, is a finite or denumerable set on which 
a size function is defined, satisfying the following conditions: 

1 . the size of an element is a non-negative integer; 

2. the number of elements of any given size is finite. 

If ^ is a class and a is an object in the class, we denote the size of a by \a\ or \a\^ when 
more than one class is under consideration. For a lattice walk w 6 'W, \w\ is the length of 
the path, or the number of steps taken. For a class Jl, we usually denote by J?l„ the set of 
objects in ^ which have size n, and by a„ the cardinality of yi„. 

Definition 8. The sequence (a„) is known as the counting sequence of the combinatorial 
class J{. 

The notion of isomorphism presents itself in this context as follows. 

Definition 9. Two combinatorial classes, J{. and S are said to be (combinatorially) isomor- 
phic, written J?l = S, if and only if their counting sequences are identical, that is, a„ = bn 
for every n > 0. This condition is equivalent to a bijection from J?l to S which preserves 
size, and we say that J?l and S are bijectively equivalent. 

It should be noted that while two classes may be equinumerous, and able to be shown 
to be isomorphic, it is not always easy or possible to find a nice or natural transformation 
between two classes that bijectively maps one to another. The enumerative information 
about a class J?I is stored in a formal power series A{t). 
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Definition 10. The ordinary generating function (OGF) of a sequence (a„)„ is the formal 
power series 

CO 

A{t) = J^a,,t". 

n=0 

The OGF for a combinatorial class Jl is the generating function for the counting sequence 
an = \^n\, with the variable t marking size. 

There are two special classes: the empty class and the atomic class. The empty class, 
denoted by fi, is the unique class (up to isomorphism) containing one object of size and 
has OGF E{t) = 1. The atomic class, denoted sometimes by Z., is the class containing one 
object of size one with OGF Z{t) = t. In some applications we differentiate between atoms: 
for instance, to represent steps in different directions. When we differentiate a NE step and 
a SE step, we would use Z.ne and Z.se- 

Just a brief note on terminology: since the counting sequence for the class ^ is also the 
coefficient sequence for the generating function A{t), we shall be using both terms when 
referring to the sequence («„). 

Note the naming convention at work here: a class ^ has subclasses Jl„ of elements of 
size n, each with cardinality a^, and generating function A{t). Also, sometimes we would 
like to find the coefficient of the t" term in the power series. We call this action coefficient 
extraction, and denote it by 

[?"]A(0 = an. 

Generating functions are manipulated algebraically as formal power series. Formally, 
they are elements of the ring of formal power series C^x^. The operations of sum and 
Cauchy product in the ring of formal power series are interpreted quite naturally as oper- 
ations on the combinatorial structures themselves. Let's investigate these operations: let 
Jl and S be a pair of combinatorial classes with no intersection. We define two binary 
operations on Jl and S: combinatorial sum (+) and combinatorial product (•). 

We define the sum of disjoint classes Jl and S like so: 

J?1 + S:=J?IUS. 

Inside our new unified class, the class of objects of size n is the union J?l„ U S„. Since 
J?l n S = the cardinality of J?l„ U S„ is simply a„ + bn- We call the new sum class C, which 
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has OGF 

C(t) = Yjittn + bn)f, 
n>0 

which is clearly the sum of the OGFs A(t) and B{t). 

Next, we define the product of ^ and S as the cartesian product: 

The elements of the product class are ordered pairs {a,/3), and the size of an ordered pair is 
defined to be the sum of the sizes of its entries. So if \a\ = 2 and \/3\ = 4, then \{a,/3)\ = 6. 
We call the new product class P. Let us look again at the subclass ^„ of objects of size n. 
The elements are ordered pairs (a,>S) with, for each k = 0, ..., n, \a\ = k and \/3\ = n - k. For 
each k there are atbn^k such ordered pairs, so the generating function for this product class 
is 

P(t) = J]\J]atbn^ 



n>Q \k=0 J 



t\ 



which is precisely the Cauchy product of the generating functions A{t) and 5(0- 

We may now extend the definition of the sum to include classes that share some ele- 
ments. Suppose that SI and S did have a non-empty intersection. Then we would simply 
take 

j?l + S:=j?luS, 

where U is a disjoint union of the two classes. We take this by representing each class as a 
product: J{. = J[x& and S = fi x S. Since & is the empty class, this transformation adds 
no weight to the class (hence the congruence asserted above is legitimate) and the disjoint 
union still corresponds to the sum of the OGFs. 

For many of our cases, we will not be finding a generating function, but a functional 
equation that the generating function satisfies. Take for example the geometric series 



W{t) = J^ Tf. 



n>0 



This is the generating function for the number of unrestricted lattice walks on the step 
set {NE, SE}. This satisfies the algebraic functional equation 

W(t) = 1 + 2tW(t), 
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which is an algebraic equation with polynomial coefficients. Thus the solution, the gener- 
ating function, is an algebraic function given by the sum of a geometric progression 



W(t) 



1 



\-2t 



Solutions to algebraic equations exist and are unique up to multiplicity of roots, and this 
fact gives us a quick way of checking an isomorphism between two classes. If two classes 
have generating functions which satisfy the same functional equation, then the coefficient 
sequences of each will satisfy the same recursion. To prove isomorphism, all that remains to 
be done is a comparison of the first few terms of each sequence. A good example where the 
isomorphism fails due to a diff'erent sequence are the neutral class, with GF E(t) = 1, and 
atomic class, with GF Z{t) = t, which both satisfy the same algebraic equation A{t)^ = A{t) 
but are clearly not equinumerous. 

Example 4.1. (Dyck paths.) We can use the theory developed so far to investigate what is 
probably the most famous example of a class of lattice paths, D, the class of Dyck paths. 
These are paths on the step set S = {NE, SE}, beginning at the origin and ending on the 
jc-axis, without leaving the first quadrant - in the nomenclature of [2] they are an example 
of a directed excursion. A Dyck path of length 16 is shown in Figure 14X1 




Figure 4.1: A Dyck path of length 16. 



These paths are well known, and the number of Dyck paths of length 2n is known to be 

the nth Catalan number, 

1 (2n\ 
C = 



n + l\ n )' 
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so the OGF for D is given by 



n>0 ^ ' 



We can use the operations of sum and product to find a functional equation satisfied by 
the OGF for D, which we may then solve and use a coefficient extraction to find the Catalan 
numbers. This is done by decomposing a Dyck path into 2 parts, using what is known as the 
inital pass decomposition. The decomposition is the following: a member of the class D is 
either the empty path or a path of non-zero length. If it is a path of non-zero length, after 
the point of contact with the x-axis formed by the beginning of the path, there will be at 
least one more point of contact. This is known as the initial pass, and is indicated in Figure 
|4.2| by a blue circle. 

Prior to the initial pass there are only two steps in contact with the x-axis, which are 
shown in blue in Figure 14.21 Between these steps is a shorter (possibly empty) Dyck path 
beginning and ending on the dashed blue line. After the initial pass there will be another 
(again, possibly empty) Dyck path beginning and ending on the x-axis. 




Figure 4.2: The initial pass decomposition of our Dyck path 



This is summarised symbolically in terms of 2), the empty class S and the atomic 
class Z as 

D = 6 + ZNE*2)*ZsE*n. 

We then use the definitions of combinatorial sum and product to interpret this in the ring of 
formal power series containing the OGF for D: 

D{t) = 1 + tD{t)tD(t) = 1 + t^D{tf. 
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Still working with D{t) as a formal algebraic object, we solve the equation, a quadratic 
in the OGF D{t), using the quadratic formula: 



-(o = i^|=^. 



We may reject the branch 



1+ Vl - 4f2 



2fi 
since when we develop its power series (by taking Newton's expansion of (1 - 4-fy^^) we 

find that the corresponding power series is singular at and contains negative coefficients. 

Thus, we find that 

1- Vl - 4/^2 

m = 3,. ■ 

which we may expand to its power series form by Newton's generalised binomial expan- 
sion. 

Finding generating functions as solutions to algebraic equations brings up the question 
of classification. As with numbers, formal power series which are the solution to algebraic 
equations are known as algebraic functions. Two examples are the formal power series W{t) 
and D(t) that we have seen already. In [|24l . Stanley gives a classification of power series 
into a hierarchy, which we shall discuss presently. 

Classification of OGFs 

Let G denote an arbitrary formal power series with coefficients taken from a field K. Then G 
is an element of the ring of formal power series K^t^. The power series G is classified 
as either rational, algebraic, transcendental D-finite or non D-finite, terms which will be 
defined in due course. In Chapter [81 we find that a family of D-finite generating functions 
belong to a smaller class known as G-series, but we don't talk about these in this chapter. 

The smallest class is that of rational power series, or those which can be represented as 
the ratio of two polynomials. 

Definition 11. A formal power series G e K^t^ is rational if there exist polynomials 

Q{t), P{t) e K[t], with Q(t) ^ 0, such that 

Pit) 
G= —-. 

Qit) 
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Example 4.2. The power series 

W(t) = 



l-2t 
is an example of a rational series. 

Note that a rearrangement of the equation in Definition [TTI gives a linear equation for 
which G is the solution. The definition can be extended to equations that are of higher 
degree in G, giving the general definition of an algebraic power series. 

Definition 12. A formal power series G e K^t^ is algebraic if there exist polynomi- 
als Po{t), Piit), ..., Pdit) e K[t], not all identically 0, such that 

Po(t) + P,it)G + ... + PAt)G' = 0. 

The minimal integer d for which such an equation holds is the degree of G. 

Example 4.3. As we've already seen, the generating function 



1- Vl - Af 

for Dyck paths satisfies the algebraic equation 

1 - D{t) + t^D{tf = 0, 
making it an algebraic series of degree 2. 

The class of algebraic series sits inside a larger class, known as dijferentiably finite 
power series, shortened to D-finite. They are also known as holonomic functions: a power 
series is holonomic if it satisfies a linear homogeneous diff"erential equation, which we shall 
see is very similar to D-finite. 

Definition 13. LetG e K^t\ If there exist polynomials /^oCO^ ■■■,Pd(t) 6 K[t] withp^it) i^ 
such that 

pAm"^'' + p.-.im'-'-'^ + ... + pi(OG' + po(OG = 0, (4.1) 

where G'^^ = ^, then we say that G is a D-finite power series. 

Algebraic functions are D-finite, by the following theorem, taken from Stanley's f24\, 
along with a sketch of the proof. 
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Theorem 4.1.1. / i2?l Theorem 6.4.6]Let u e ^[[x]] be algebraic of degree d. Then u is 
D-finite. 

Proof. Let u be as stated. Since u is algebraic, there is some degree d polynomial Pix,y) e 
K[x, y] such that P{x, u) = 0. Thus 



J dP{x,y) 

= -—P(x, u) - 



dx dx 

Since P{x, y) is of minimal degree, 

dP{x,y) 



,dP{x,y) 
+ u 



y=u 



dy 



y=u 



dy 



^0, 

y=u 



and so 



u = 



dP{]c,y) I 



dP(x,y) 



dy 



y=u 

That is, u' is a rational function in x and a, so is an element of the field K{x,u). We 
can continue inductively to show that for any k > 0, u^'^^ 6 K{x, u). Now, u is algebraic of 
degree d, so AimK^x) K{x, u) = d, and therefore u,u', ..., u^''^ are linearly dependent over K(x), 
giving an equation of the same form as Equation (14.11 ). proving that u is D-finite. n 

Let's consider a brief example of a D-finite function that is not algebraic. 

Example 4.4. An example of a D-finite function is G(0 = fe\ since G'(0 = (l + j)G. How- 
ever, f-e' is not algebraic. This shows that the algebraic functions form a proper subclass of 
the D-finite functions. 

Finally, if a power series G e K^t^ is not D-finite, we simply say that G is non D-finite. 
Proof that a generating function is non-D-finite can be difficult, see [|22| for two lattice path 
examples from a family for which classification is ongoing. 

The two expressions 

W{t) = ^ 



D(t) = 



\-2f 

1- Vl-4r2 

If 
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raise an interesting point about these generating functions: we may view them as standard 
complex analytic objects by assigning complex values to the variable t. Both series will 
converge to analytic functions for \t\ < ^. The exploitation of the complex analytic proper- 
ties of the generating functions for enumerative results will form the focus of Section 14.31 

The discussion so far has focused on ordinary generating functions, those with a single 
variable. Lattice paths require something more flexible than this, and Section |4]2]introduces 
the required ideas. 

4.2 Multivariate generating functions 

The formal power series framework is easily adapted to a multivariate model, where mul- 
tiple parameters are introduced. When working on the main problem, we use a trivariate 
power series with parameters keeping track of the end point of a walk. We show a simpli- 
fied version: a bivariate generating function (BGF) which marks both the length of a walk 
and the y-coordinate of its endpoint, which we call the stopping height. 

Example 4.5. We look at the class "W of walks taken on the step set {NE, SE}, but unre- 
stricted in where it may visit and end. Walks can be built recursively from shorter ones by 
simply appending a new step, and this can be done in two ways. Symbolically, we have 

nV = £ + nV • Z;v£ + ^ • ZsE, 

which we may give as the familiar functional equation 

Wit) = 1 + 2tWit). 

We keep track of the stopping height, and do so by introducing a parameter, u. Note that 
to keep track of height, we must increment the exponent of w by 1 for an upward step, and 
decrement by 1 for a downward step. To capture this, we move straight into the ring of 
formal power series, rather than using a symbolic specification. The only diff"erence is that 
now K = C(m) because the walks are allowed to visit points with a negative y-value. The 
corresponding algebraic equation is 

W*it, u) = l+ utW*it, u) + - W*{t, u), 

u 
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which we may solve to find the MGF as the geometric series 

1 



W*(t,u) 



i-t(u+iy 



where t marks length and u marks stopping height. We can then perform a coefficient 
extraction with respect to t to find a Laurent polynomial in u 

[nw*(t,u) = iu + -] , 

where the coefficient of u'' is the number of walks of length n with stopping height k. To 
perform the coefficient extraction, we use the following fact: 

[u'-']f(u) = [u']u'f{u). (4.2) 

This allows us to simplify things a little when expanding the binomial. Note that the range 
of values possible for k is the set {-n, -n + I, ..., n - l,n}, but to simplify calculations we 



multiply by u" (as in Equation 14. 21) so the range of interest in the new function is [[0, 2n]]. 
So, for k e [[0, 2n]] we have 



This is a binomial with a quadratic term, so for ^ = 1 (mod 2), we have 

[u'']iu^ + 1)" = 0, 
and fork = (mod 2), we have 

So for k e {0, 1, 2, ..., 2n} the coefficients for our MGF are 
[t"u^-"]W*{t,u) = 



Mr. .2 , 1^n_/" 

k 

2 



, I 0, ifyt= 1 (mod 2) 



(Jj ifk^O (mod 2). 



The classification of formal power series also applies in the multivariate case. Now we 
take G 6 K[[ti, ..., ?„]], and classify it at as rational, algebraic, D-finite or non D-finite over 
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the field K. The first definition is the same as the single variable case: G is rational if there 
exist polynomials P,Q& K[ti, ..., tn\ with 2^0 such that QG = P. 

The definition for algebraic is also largely the same. If there exist polynomials Pq, ...,Pd e 
K{t\, ..., tm] not all such that Pq + PiG + ... + PdG'^ = 0, then we say that G is algebraic. 
If d is the smallest integer such that this is true, then G is algebraic of degree d. 

Finally, the generalisation of the D-finite definition requires more. A power series G e 
K[[t\, ...,?„]] is D-finite over K(t\, ...,?„) if for 1 < i < m,G satisfies a system of non-trivial 
partial differential equations of the form 

where Pjj e K[ti , ...,tm\ and J, is the order of the partial differential equation in ti. 

4.3 Coefficient asymptotics 

For classes for which an exact counting formula is known, real analysis is usually enough 
to obtain asymptotic expressions for the counting sequence, as we show with the Catalan 
numbers. 

Example 4.6. We show how to get asymptotic expressions from a formula by applying 
Stirling's approximation to the formula for the Catalan numbers. Recall from real analysis 
that Stirling's approximation says that 

n\ ~ \2Tm\-\ . 

We apply this to the formula 

1 (In 
C 



n + \\n I 
while noting that n + l~n.So, we get 



n + l\n I 

1 (2nf"e-^" ^[Aml 

02k 



Inn^ 
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Next is the class of problems for which a generating function in a closed form is obtain- 
able, but no closed form expression is available for the counting sequence. It usually suffices 
to have an expression of the generating function as an analytic function, not a formal power 
series, to obtain asymptotic expressions. The analysis depends on the local properties of the 
generating function at its dominant singularities. The exposition of this theory is the focus 
of this section. To continue, we recall the definition of asymptotic growth. 

Definition 14. Let a : N ^ N be a function (such as the coefficient sequence of a generating 
function associated with a combinatorial class). We say that a„ grows like or has asymptotic 
growth g{n), ^ : R"^ ^ R"^, denoted 

a{n) ~ g{n) 

if for every < e < 1 we can find a sufficiently large A^ such that forn>N the following 

holds 

a{n) 

1 - 6 < — — - < 1 + 6. 

g(n) 

As we can see in Example 14.61 and for our rational function example, the coefficients /„ 
of a generating function F belong to a general asymptotic regime. That is 

[t"]F(t) - /S^ein), 

where /3" is an exponential growth factor modulated by a subexponential factor 0{n). For 
Dyck paths, the exponential growth factor is 4, but recall that we count only even length 
paths so we may replace 4" by 2^", and for the rational function, the coefficients are 2". 
Recall that in both cases, the radius of convergence is ^. This (tenuously) suggests a link 
between the radius of convergence of a generating function and the exponential growth of 
the coefficient sequence, and this tums out to be the case. 

We begin by quoting from [ | ir | the two principles of coefficient asymptotics. 



The first principle of coefficient asymptotics: The location of a func- 
tion's singularities dictates the exponential growth of its coefficients. 
The second principle of coefficient asymptotics: The nature of a func- 
tion's singularities determines the associated subexponential factor. 
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So, rather than viewing the generating functions as formal algebraic structures, we begin 
to view them as analytic functions of a complex variable. These can in turn be viewed in 
two equivalent forms: as convergent power series or as differentiable functions. The first 
form directly relates to the use of generating functions for enumeration. 

We assume a basic understanding of complex analysis: the definitions of analytic, holo- 
morphic and meromorphic functions are left to more focused texts [IT9|, as well as the 
definition of a singularity. In order to use information from the singularities in the asymp- 
totic expressions for the coefficients of the generating function, we need to relate the global 
properties of the function near - the asymptotic expansions - to the local properties of the 
function away from zero - the singularites. 

Cauchy's residue theorem and Cauchy's coeflicient formula are both examples of the- 
orems which relate global and local properties of functions through the use of contour 
integrals, with contours taken around the singularities of a function. In practice, Cauchy's 
coefficient formula is often used to find bounds on the asymptotic growth of a generating 
functions coefficient sequence, as in [|2l Theorem 3]. Given that we are interested in finding 
exponential growth factors, we need only consider singularities of a generating function 
which lie at the boundary of the disc of convergence. These are called dominant singular- 
ities. The following theorem of Pringsheim significantly simplifies the search for these in 
the case of generating functions. We omit the proof, but it can be found in [fTTI . 

Theorem 4.3.1. (Pringsheim' s theorem). If f{t) is representable at the origin by a con- 
vergent series expansion that has non-negative coefficients and a radius of convergence R, 
then the point t = Ris a singularity of f{t). 

This tells us that a dominant singularity of a generating function will lie on the real 
line, and allows us to only consider the real singularities of a generating function when 
determining the radius of convergence. In order to show the relation between the radius of 
convergence (or location of the dominant real singularity) and the exponential growth of a 
counting sequence, we first need a definition. 

Definition 15. We say that a number sequence (a„) is of exponential order K" which is 
denoted as 

a„xK" 
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if and only if 

limsup|a„|" = K. 

n— ►oo 

The relation a,, x K" expresses both an upper and lower bound: for any 6 > we have: 

1- M«l >i o (^ ~ ^^"' ^^^^ ^^' l^nl exceeds (K - e)" infinitely often; 

2. \cin\ <a.e {K + €)", that is, \a„\ is dominated by (K + ef almost everywhere. 

Often, our sequences have a subexponential growth factor d(n) as well. Note that such 
a factor would be dominated almost everywhere by an increasing exponential (1 +6)", and 
bounded from below infinitely often by any decaying exponential (1-6)". Thus 

limsup|0(n)|" = 1, 

n— »oo 

and if a„ = K"6(n), it still satisfies a„ x K". We are now ready to prove the exponential 
growth formula. 

Theorem 4.3.2. tUl, Theorem IV.7]( Exponential growth formula). If f(z) is a generat- 
ing function, analytic at 0, with a power series representation involving only non-negative 
coefficients, and R is the modulus of the dominant real singularity, that is 

R = sup{r >0\fis analytic at all points ofO < t < r}, 

then the coefficient f„ = [t"]f(t) satisfies 

Proof. Let R be as in the statement of the theorem. Pringsheim's theorem tells us that R is 
the radius of convergence. Now, we must prove the exponential growth relation. Since R is 
the radius of convergence, the power series will converge at any point just to the left of R, 
that is, for 7? > e > 

f„{R - e)" ^ 0. 

Li other words 

f„(R -€)"<! 



CHAPTER 4. ANALYTIC COMBINATORICS 31 

for all n > N, given N sufficiently large. This gives 

\f/^<(R-er' 

almost every where. For property 2, let e > 0. Then/„(7?+6)" cannot be a bounded sequence, 
otherwise the power series 

J]UR + e/2r 

n>Q 

would converge. Indeed, suppose that \fn(R + e)"| < B for some bound B. Then 

\f„\(R + e/ir = \f„\{R + eT i^^^] ^ 0, 

since \f„\(R + e)" is bounded and 

R + e I 
Therefore /„''" > (R + t) ' infinitely often, giving the relation 

n 

This proves the first principle of coefficient asymptotics. The second principle, relat- 
ing the subexponential growth of the counting sequence to the nature of the singularity, is 
dependent on the Newton expansion of the generating function near the singularity. We 
illustrate the previous theorem with some examples. 

Example 4.7. The first example is of lattice paths with the same set of steps as Dyck paths, 
but no longer restricted to lie in the first quadrant: they can now joumey as far south as 
they would like and end at any place they can reach. We'll call the class of such paths S. It 
should be clear that since there are no restrictions, the generating function is 



B(t) = J]2"f = j^, 



n>0 



This clearly has a simple pole att = ^, and no other singularities. Applying the exponential 
growth formula, we find 

bn = reBin), 



CHAPTER 4. ANALYTIC COMBINATORICS 32 

where Osin) is a subexponential growth factor. Note that in this case 0B(n) = 1, since the 
counting sequence is exactly 2". 

The second class is the class of Dyck paths, D. We already know the generating func- 
tion: 

1 (2n\ 2„ 1 - Vr^4?2 






There is an apparent pole at 0, but development of the generating function as a series will 
show that this is removable. The dominant singularity in fact lies at ? = ^, as in the previous 
case. This is the branch point of the square root portion of the function: the function Vz 
cannot be unambiguously defined in a neighbourhood of Z = 0. Applying the exponential 
growth formula gives 

dm = 2^"eD{2n). 

This time, due to the nature of the singularity the subexponential growth factor is a non- 
trivial function of n, which reflects the restriction of the paths to lie in the first quadrant. 
We already applied Stirling's approximation to the formula for Catalan numbers in Example 
so we know the asymptotic growth of the function 9 to be 

1 



0Din) 



Inri^ 



The final class we consider is the class of Motzkin paths A\. Like Dyck paths, these are 
restricted to lie in the first quadrant and begin and end on the jc-axis, but use the set of steps 
shown in Figure 14.31 



^ 



Figure 4.3: The steps for Motzkin paths 

Like the Dyck paths, we can decompose Motzkin paths. Again, we either have an empty 
walk, or we take a step. The first step is either horizontal or an upward step. If the first step 
is horizontal, what follows will be a Motzkin path. If the first step is upwards, we use 
the first pass decomposition as before. This gives the functional equation satisfied by the 
generating function 

M{t) = 1 + tM(t) + t^M{tf. 
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We solve this algebraically to find the algebraic function the generating series converges to 
inside its radius of convergence: 



M(,) = i-'^^i^jm-m 

Since the generating function is analytic at zero, we know that the apparent pole at ? = is 
removable. This makes the dominant singularity ? = |, giving 

M„ = rOMin). 

Again, this time Om is non-constant, and o(l). An application of Stirling's approximation 
to one of the manifold Motzkin number formulas 






feO 



where Q is the ^th Catalan number, shows us that 

3V3 



dM(n) 



The second principle of coefficient asymptotics states that we must look to the nature 
of the singularities for the subexponential growth of a counting sequence. For many of the 
cases we consider, a closed form expression for the function is as yet unknown. This makes 
an application of Stirling's approximation as in Example 14.61 impossible. Therefore, we 
must turn to the classification of the generating function for asymptotic growth templates. 
For general algebraic series. Wimp and Zeilberger give us a method in [|29l for finding an 
asymptotic growth template for the coefficients: 

an ~ K/3"n" 

where /5 = ^ and k are algebraic over Q and a e Q \ {-1, -2, ...}. That is, a constant, 
an exponential part and a polynomial part. The more general class of D-finite generating 
functions also has a template that we may find through the methods of f29^. This one is 
more complicated: 

an ~ /<}e'V^("*"")n''(log(^))^ 

with l,r eN, K,a algebraic over Q and P a polynomial. 
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As shown in [4J, all the D-finite series we consider are examples of G-series, so as 
mentioned in the introduction we actually have the simpler growth template 

an - /c-yS"n"(log(n))^. 

The only difference between this and the algebraic case is the logarithmic factor, the expo- 
nent of which, 7, is found to be zero in all D-finite cases that we consider (through methods 
given in both [|29l[Tni ). Thus, the growth template for models of interest is the same as the 
algebraic case. 

This is all we say on the subexponential factors, since the exponential factors /3 are 
what we are more interested in. The reasoning behind this was given in Chapter |3] the 
exponential growth factors are closely linked to the entropy of the physical system being 
modeled. This concludes the background information for the thesis. The next section will 
survey four different lattice path models, beginning with the most simple to analyse in 
Chapter m 



Part II 



Planar lattice paths with small steps 
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Chapter 5 



Unrestricted walks 



We begin the second part of this thesis by considering the most general case of planar 
lattice paths with small steps, those completely unrestricted in where they visit. One can 
quite trivially enumerate this model for any step set, but we wish to pedagogically expose 
some of the developed machinery at work on something within the grasp of any reader. 

Section [5TT] will introduce the class of unrestricted planar lattice paths, giving the ordi- 
nary and multivariate generating functions and producing the functional equation satisfied 
by them from a symbolic specification. The functional equation is quite easy to solve, and 
makes classification in Section 15^ simple. Finally, we exactly and asymptotically enumer- 
ate the class in Section [531 

5.1 Definition 

Let S c S be a set of small steps, and consider lattice paths with steps from S, beginning 
at the origin and unrestricted in the regions of the plane that they may visit. This is the 
class 'Ws of unrestricted lattice paths on S. A typical example was shown in Figure IZ2l 
Taking w(n) to be the number of walks in 'Wz of length n, we define the ordinary length 
generating function as 

We can also keep track of the endpoint of a walk. If we let w{i, j; n) be the number of 
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walks of length n ending at the point (/, j), then we can define the MGF for "Wq as 

Ws(x,y;t)= Yj w{ij\n)xyf. 

(',;)eZ2 
n>0 

Here, x and y mark the respective coordinates of the endpoint. Note that W{t) = V7s(l, 1; 0- 
Since this class is completely unrestricted, the recursive structure of the walks is obvi- 
ous: all walks of length n + \ can be created by appending a step from S to the end of a walk 
of length n. We can translate this into a specification, using the operations of combinatorial 
sum and product: 

nVs = fi + S • Z • ^s, 

where the summation is combinatorial in nature. In the next section, we will translate this 
into the ring of formal power series and classify the generating function. In this ring, the 
product of the step set and a class will translate into the product of the inventory and the 
generating function for that class. 

5.2 Classification 

The MGF (and hence the OGF) for 'Wq is easily shown to be a rational power series. Using 
the techniques from Chapter |4l the specification given for "Wg in Section 15.11 is translated 
into a multivariate functional equation, translating a combinatorial sum and product into 
the operations of sum and Cauchy product in the ring of formal power series. This gives 

Ws(x, y\t) = \+S{x,y)-t- W^ix, y; t), 

which we may solve algebraically to find the rational function 

1 



WQ(x,y;t) 



l-Six,y)t 

which the MGF will converge to inside its radius of convergence. Setting x = y = I gives 

1 



W^it) = 



i-mt 

as the rational function for the OGF. 
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Example 5.1. Let S = A. Then S(x,y) = y + — + -, and the functional equation satisfied 



by W(5(x,y;t) is 



t tx\ 
WQ{x,y\ t) = l+\ty + — + —\ WQix,y\ t), 

xy y 



which gives the rational MGF 

We(x,y;t) 



and rational OGF 



l-3t 

5.3 Enumeration 

Note that we can count these directly using the same observation that gave us the recurrence 
in Section |5l2l at each step we have |S| choices for the next step, so the number of walks is 

w(n) = I SI". 

This allows us to define the OGF immediately as the power series 

Ws(t) = J]\ert\ 

n>0 

which is a geometric series, so has sum 



i-m 



fort < 7^. Incidentally, this is the radius of convergence of the power series Ws(0, and Sec- 
tion |4.3l gives us a method for finding the asymptotic growth of the sequence w{n) using this 
radius of convergence. But since rational generating functions have geometric power series 
and thus their coefficient sequences grow purely exponentially, the asymptotic expression 
is the same as the exact expression. 



Chapter 6 



Directed walks 



Directed paths are those in which all steps have a privileged direction of increase. In this 
instance, we take all steps to be x-positive. Since this direction of increase forces self 
avoidance, this makes these paths naturally suitable for modeling of linear polymers subject 
to a force, as discussed in Chapter|3] We primarily use them in our analysis of quarter plane 
walks in Chapter [8] 

In this chapter, we discuss results from Banderier and Flajolet's paper [|2l on directed 
paths which stay above the x-axis, known as directed meanders. The introduction of a 
boundary affects the enumerative results for walks on some step sets. This is due to a 
parameter known as the drift of the walk, which we discuss in Section 16. 3[ A common 
example of a directed meander, known as a Dyck prefix, is shown in Figure 16711 



> 



Figure 6.1: A Dyck prefix of length 10 on S = {NE, SE}. 
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First, a brief note about the step sets used in this chapter. In L2J, walks are not restricted 
to just small steps. Steps may in fact have any j-coordinate as long as they step one unit to 
the right. Since the quarter plane models considered in Chapter [8] are all on sets of small 
steps, we make the restriction here also. As in the other models, restricting to small steps 
significantly simplifies the computations necessary (for example, compare Theorem 16.2.11 
here and Theorem 2 in [2J). Future work on these problems include a generalisation of the 
results to larger steps, see Chapter |9]for a brief discussion. 

6.1 Definition 

We call the class of directed meanders 'F, as in |l2l. In order to modify our definition from 
Chapter |2l we add the restrictions that each step has an .x;-value of 1 and that each point 
visited by a path lies in the half plane y > 0. For walks with small steps, this shrinks the 



largest step set to that shown in Figure 16.21 



Figure 6.2: The set of directed, small steps 

For a given set of directed small steps S, the inventory S (x, y) contains superfluous 
information: we know that the x value of every step is 1 . Thus we may omit the x-coordinate 
information. In light of this, we define a new directed inventory 

P(y):=S{\,y). 

Example 6.1. Let S be the full set of directed small steps shown in Figure |6^ Then S has 
the inventory S(x,y) = xy + x + x/y and directed inventory P(y) = y + I + l/y. 

For a given set of directed small steps S, let f(j; n) be the number of directed meanders 
with endpoints having >'-coordinate, or final altitude, j and length n. Then the MGF for 
directed meanders on S is 

F^{y;t)=Y,fU\n)yh\ 
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Again, these walks can be built recursively: we create directed meanders of length n + I 
by appending steps from S in all possible ways to meanders of length n. The phrase 'in all 
possible ways' is a little vague: what we mean by this is that if a walk of length n has its 
endpoint on the x-axis, we may not create a walk of length n + 1 by appending a j-negative 
step. This last point is difficult to reflect symbolically in terms of combinatorial sum and 
product, so we immediately define a functional equation based on the recursion: 

Fdr, 0=1+ tP(y)Fdr, t) - {j^"} [P(y)] ■ tFs{0; t). (6.1) 

Note that the final term represents removing the 'impossible ways' to increase the length of 
a shorter walk. Rearrangement of the functional equation gives the expression 

1 - {y<"}?P(y)Fs(0; t) 

^^^^'^^ = TlPOO • ^^-'^ 

Example 6.2. Take S = {NE, SE}, the steps of a Dyck prefix. Then the directed inventory 

of S is 

P{y)=y+-, 

y 

and the functional equation satisfied by F^{y; t), the BGF for the class of Dyck prefixes, is 

F^iy; t)=\+t\y + -\ F^iy, t) - -Fs(0; 0- 

As in Equation (16.21) , we may rearrange the functional equation to find the expression 

l-iFdO;t) 
Fs(y;t) = 



1-K^ + i) 

The expression (16.21) is for the MGF of directed meanders on a step set S in terms of 
the OGF for meanders ending on the x-axis. Unless we know the classification of Fs(0; t), 
we are unable to immediately classify F^iy, t) as we could for the whole plane model. 
In Section [6^ we shall classify the generating functions for directed meanders with small 
steps by solving the functional equation explicitly through the kernel method. 

6.2 Classification 

We continue with expression (16.21) 

^ ^ ^^ l-{y<''}tP{y)F^{S)-t) 
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for the MGF of directed meanders. The quantity 1 - tP{y) found in the denominator of the 
above expression is important enough to warrant a name. 

Definition 16. Let Ts be a class of directed meanders taken on a set of directed small 
steps S with directed inventory P(y). The quantity 

1 - tP{y) 

found by rearrangement of Equation (|6.11 ) is called the kernel. We define the kernel equation 
to be 

y - tyP(y) = 0, 

where the left hand side of the kernel equation is the kernel made entire. 

We work with directed small steps so the kernel equation is at most quadratic in y. Solv- 
ing by the quadratic formula will yield two solutions, both of which are algebraic functions 
of t. When developed as power series, one shall be convergent at zero and have only non- 
negative integer coefficients and the other will be singular at 0. Denote the solution with 
non-negative integer coefficients as yo{t). This function plays a significant role in the solu- 



tion of Equation (16.11) . 

Example 6.3. We continue Example 16. 2 [ which showed us that the functional equation is 

Fs(y; t) = l+t(y+ l/y)Fdy; t) - t/yF^iO; t), 

and rearrangment will give 

\-t/yFdO;t) 

l-t(y+ l/y) 

We make the kemel entire, and set it equal to zero to get the kernel equation 

y - ty^ - t = 0, 

which we may solve using the quadratic formula: 

_ -1 + VI - 4(-0(-0 _ 1 ± Vl - 4f 
^ ~ 2^0 ~ 2t ■ 

Development of Vl - 4t^ will show that 

1+ Vl - 4f 

y^'^ = — 2t — 
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is not convergent at t = 0, so 

^^ 1- Vr^4f2 

'^^'^ = 2t 

is the desired solution to the kernel equation. 

We now move on to a proof that for all directed meanders taken on small steps, the 
generating function is algebraic. In the process, we obtain a formula for Fs(y; t) in terms 
of the kernel and the power series yo{t)- 

Theorem 6.2.1. ^ Theorem 2] For a set of directed, small steps £, the BGF of meanders 
(with t marking size and y marking final altitude) relative to S is algebraic. It is given in 
terms of the function yoit), the solution to the kernel equation, by 

The counting generating function is found by taking y = 1, giving 

Proof. The fundamental equation in the form of Equation (|6.11 ) looks underdetermined. It 
involves two unknowns: the bivariate F^iy, t) and the univariate Fs(0; t). The main idea of 
the kernel method is to bind t and y in such a way that the left hand side vanishes. 
First, we remove the negative exponents of y from Equation (16.11 ): 

yF^iy; t)=y + tyP(y)F^(y; t)-y ({3^<°}/'(y))Fs(0; t). 

Since Piy) has a negative degree of 1 , {{y'^^}P{y)) is simply a Laurent monomial in y. Denote 
it by riy). We remove the pole at joCO by multiplying both sides by the kernel: 

FeCj; t)(y - tyPiy)) =y- tyr(y)Fe(0; t). 

We want to force the left hand side to vanish. By restricting ? to a small enough neighbour- 
hood of the origin contained in |;^| < 1/P(1), we can ensure that \yo{t)\ < 1 (recall 3^0 has 
non-negative coefficients), ensuring that F^iy; t) converges upon substitution. 

Thus, the left hand side of the equation vanishes providing us with an equation in a 
single unknown, the function Fe(0; t): 

= yo{t) - tyo(t) ■ r(yo)FeiO; t). 
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Thus Fz(0; t) is an algebraic function expressible rationally in terms of the algebraic func- 
tion joCO- 

Next, we make use of an observation of Mireille Bousquet-Melou [HI. Define 

N{t,y):=y-tyriy)F^{0-t). 

This is by the substitution in the previous paragraph a linear polynomial in y with its root 
at the small branch, yo{t). Thus, we can simply write it as 

N{t,y)=y-yG{t). 

Thus, the constant term is simultaneously -jo(0 and -tyr{y)F q{0\ t), which we may solve 
for Fs(0; t). Thus, N{t,y) is fully determined in terms of the small branch, and we get the 
MGF for meanders from the fundamental functional equation made entire: 

F ( ■t)= ^^^'^^ = y-yo^^^ 

This gives the OGF of meanders 

^ 1 - tP(\) 
Since the algebraic functions form a field, and N{t, y), P{y) and yo{t) are algebraic functions, 
both the MGF and OGF for directed meanders are algebraic. n 

Example 6.4. We may apply Theorem 16.2. II to Example 16. 3 1 to produce a generating func- 
tion for the class. Recall that the directed inventory for the step set is P(y) = y + - and the 
solution to the kernel equation is 



^^ 1- Vl-4f2 

^°^^^ = 2t • 

Theorem 16 . 2 . 1 1 tells us that the bivariate generating function for the class of Dyck prefixes 

is 

2ty-H- Vr^ 

and the OGF is 

2t-\+ Vl - 4-f 

^-<" = (l-2,)2, ■ 
This concludes the classification of the generating functions associated with directed 
meanders. In the next section, we turn to the enumeration of these objects, both exact and 
asymptotic. 
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6.3 Enumeration 

Now that we have a method for finding the generating functions associated with directed 
meanders, enumeration seems an easy task. For closed coefficient formulae, we can simply 
apply the generalised binomial theorem and other coefficient extraction techniques to the 
generating functions given by Theorem |6.2.1| We would also like a corresponding theorem 
about the asymptotic behaviour of the coefficient sequences. As is well understood, the 
asymptotic behaviour of counts is closely related to the singular properties of the corre- 
sponding generating functions [ITTI . The application of the kernel method in Theorem 16.2. II 
gives a factorisation of the generating function under which the singular forms become 
manageable. This uses a small amount of analysis in the vicinity of what is known as the 
structural radius, p. 

Lemma 6.3.1. /El Lemma 2] Let P{u) be the polynomial inventory of a simple step set S. 
Then there exists a unique number r, called the structural constant, such that 

P'ir) = 0, T > 0. 

The structural radius is by definition the quantity 

1 



p:-- 



Pir) 

Proof. Differentiating P twice, notice that P"{x) > for all x > 0, so viewed as a real 
function, P(x) is strictly convex. Since it satisfies P{0) = P{+oo) = +00, there must be a 
unique positive minimum attained at some t, thus /"(r) = 0. n 

Finally, we can get to the asymptotic results for directed meanders. There are three 
cases, distinguished by a quantity known as the drift. 

Definition 17. []2l Definition 5] Given a simple walk with polynomial inventory P(y), the 
drift is the quantity 

6 = P'il). 

In the combinatorial case, with or without weighted steps (see Chapter |2] for how to 
weight a step set), the drift is then the sum of all the y-values of the steps multiplied by 
the weight of each step, which is an indicator of the 'tendency' of the walk to increase or 
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decrease in altitude. For a probabilistic case, the drift is the expected movement in the y 
direction of any single step. Note that for a symmetric step set, the drift is 5 = and the 
structural constant is r = 1. 

Theorem 6.3.2. [2^ Theorem 4] Let Q be a simple set of small steps, with polynomial 
inventory P. The asymptotic number of meanders taken on S depends on the sign of the 
drifts = /"(I) as follows. 



6 = 0: [f^WsiUt) - vo^(l + ^ + | + ---) 

^0 ■- -^^P"(i)' 



6<0: [t"]Fdht) - y-I^{l + U'j + ...) 

y± .= _ hf^ I • 

^0 • y p"{t) p(t)-p(1)' 



6>0: [?"]Fs(l;0 - ^o/'d)" + v^^ (l + J + f + • • •) 

^o:=i-jo(Pi), Pi = par'. 

We can view these formulae intuitively. For the positive drift cases, a fraction close 
to ^0 of the unconstrained walks is a meander, which makes sense since the walks have 
a natural tendency to gain in altitude. For negative drift, most paths tend to head down, 
and so the fraction of unconstrained walks which are meanders is exponentially small, 
roughly {P{r)/P{\))". For zero drift, the proportion of walks which are meanders is as 
large as 1/ V^, with the walks tending to oscillate near the horizontal axis. 

Proof. The discussion is based on an analysis of the generating function F^(\;t), in its 
factorisation 

Fs(l t) = 

provided by Theorem |6.2.1| The method of proof examines the locations of the singularities 
and zeroes of the numerator in relation to the zero 1 /P(l) of the denominator. 

In the case that (5 = 0, we have P'(l) = 0, r = 1 and p = l/P{l). We must do some 
analysis of the solution to the kernel equation yo{t) in the vicinity of p. First, the kernel 

equation can be put into the form 

1 
t = 



Piy) 
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Since the inverse of an analytic function at a point where the derivative is nonzero is ana- 
lytic, the above relation is invertible analytically in the neighbourhood of any point t such 
that P'{t) 4^ 0. The contrapositive then must hold: a singularity will occur at any value ^ 
such that /"(^) = 0. 

When 5 = 0, /"(I) = by construction while /"'(l) > 0. Thus the local form of 
t = 1 /P(y) reads 

t=^^-^P"{l)(y-lf + 0{(y-in 

This is readily inverted, leading to two local solutions 



Therefore, 1 -yoit) contributes a term of the form (1 -t/pY^^ att = p, while the denominator 
(l-tP(l)) has a simple zero there. Globally, the singularity of Fe( 1 ; f) is of the square root 
type, giving the result. 

For a negative drift meander, we have P'{1) < 0, so we must have r > 1 since P'{u) 
increases from -oo to +oo when y ranges from 0^ to +oo. Withp = 1/P(t) (the structural 
radius) and pi = 1 //'(I), we have pi < p. In this case, the prefactor (1 - tP{l))~^ has a pole 
at pi. This pole is cancelled, however, by a zero in the numerator induced by (1 - yo(t)) 
(since yoipi) = 1) making pi a removable singularity of Fe(l; ?)• This puts the dominant 
singularity of Fs(l; at p, and it is of the square root type. 

Finally, for positive drift, r < 1, so that the prefactor induces a pole at pi = !//'(!) 
before 1 - yo{t) becomes singular. The argument conlcudes by subtracting singularities, 
since the function 

Fs(i;t)-- -— , pi : = 



l-zP(l)' "^^ ■ P(l)' 
now has a dominant singularity of the square root type at p. n 

Example 6.5. We complete this chapter's running example, obtaining asymptotic results. 
Dyck prefixes have a zero drift step set, since 

P'(y) = 1 - ^ ^ /"(I) = 0. 

In order to apply the results of Theorem 16.3.21 we need the second derivative of the inven- 
tory: 

P"(y) = -r 
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So, by the aforementioned theorem, 



This completes the material on directed lattice paths with small steps. The following 
chapter returns to undirected paths, starting with walks confined to a half plane. By showing 
a bijection between directed walks and half plane walks, we set up a useful enumerative 
result for our analysis of quarter plane walks in Chapter |8] In Chapter |9l we discuss how 
we can use results on directed walks with more general step sets for a generalised version 
of the results given in this document. 



Chapter 7 

Walks confined to the half-plane 



We turn back to undirected walks with small steps, but begin restricting the region. By 
adding a boundary to the unrestricted case considered in Chapter|5]we can begin to see how 
these combinatorial structures could model more than just polymers. For instance, a single 
boundary model could represent the interaction of a particle with some kind of container. 
While in general any line in R^ could be used as the boundary of a half plane, we will only 
discuss the case where the jc-axis is used. 

7.1 Definition 

Fix a set of small steps S and take the class 'Wq of planar walks on S. We make the 
restriction that all elements of the sequence of points (po, ..., p„) for each walk of length n 
must lie inside the half plane y > 0. We call this subclass of walks 'Hs and seek h(n), the 
cardinality of 'Hn. Shown in Figure ITTl is a typical example of one of these. 

Given a step set S with inventory S{x,y), we can use recursion to find a functional 
equation that the generating function Hs(x,y;t) satisfies. It is similar to before: a half 
plane walk is either the empty walk or a walk of shorter length incremented by adding 
a new step. We must take care here, however, and be sure that we don't wander out of 
the allowed region. We can exculde this possibility by considering the generating function 
Hsix, 0; t), which counts walks ending on the x axis. 



49 



CHAPTER 7. WALKS CONFINED TO THE HALF-PLANE 



50 




11,S = X 



Figure 7. 1 : A walk of length 1 1 from "Ks 

So, the functional equation which H(^(x,y; t) satisfies is 

H^{x,y- = 1+ tS{x,y)Hdx,y; - {y^yS{x,y)H^{x,Q- 1). 

Notice the similarities to the functional equation for directed meanders? This is because 
these two classes are equinumerous when we make x = \, referred to in Chapter |6]as taking 
the directed inventory, given as 

S{x,y)^P{y) = S{\,y). 

This brings us to a new definition. 

Definition 18. Fix a step set S, with inventory S {x, y). When we take the directed inventory 
of S, P{y) = S{l,y), we obtain a new, directed step set Sp which is in general weighted 
(has multiple steps in the same direction), which we call the horizontal projection of S. 

We illustrate this new definition and the bijection it gives in the following example. 

Example 7.1. We show in Figure W72\ the weighted step set found by taking the directed 
inventory of S = {NE, SE, SE, NW}. This is Sp, the horizontal projection of S. 

Now in the directed model, every jc-negative step is a red step and every x-positive step 
a blue one. This gives the bijection between 'Hq on S and T(Zp on Sp, and we show an 
example of the correspondence in Figure 17.31 

This significantly simplifies classification and enumeration for us, as will be seen in the 
following two sections. 
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Figure 7.2: Taking the directed inventory of S leads to the horizontal projection Sp 
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Figure 7.3: A walk from "Hg and its horizontal projection in f^p 

7.2 Classification 

We begin the section on classification by showing the isomorphism between "He relative to 
a step set S with inventory S {x, y) and T relative to the horizontally projected step set Qp 
with directed inventory P(y) = S{l,y). 

Proposition 7.2.1. The class Hq of undirected half plane walks on a step set S and the 
class Tqp of directed meanders on Qp, the horizontal projection ofQ, are in bijection. 

Proof. Fix a set of steps S, with inventory S{x,y). We know that the generating func- 
tion HQ(x,y; t) for the class "Ke of half plane walks taken on S satisfies 

Hq{x, y\t) = l+tS {X, y)H^{x, y; t) - {y^^}tS {x, y)Hs(x, 0; t). 

By taking the horizontal projection, we are setting x = I, giving the new functional equation 

//e(l,y; 0=1+ tS(l,y)Hzi\,y; t) - {y^''}tS{l,y)H^{l,Q- 1). 

The inventory is now a function of y only, and can be represented by the inventory of the 
horizontal projection, S^. Making the replacement P{y) = 5(l,y), we find 



//s(l,3;;0 = \+tPiy)H~^i\,y-t)-{y<^}tP{y)H^i\,0-t). 
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Now, let Fzpiy, t) be the generating function for directed meanders on the set of steps 
P, the horizontal projection of S. Recall from Chapter |6] that F^^iy; t) satisfies 

F^^iy; = 1+ tP(y)Fs,(y; t) - {y^''\tP{y)F ^^^^ t), 

which is simply the functional equation for the walks in 'H^ taken on the original step set 
S with the x information suppressed (or with x set to 1). Thus, HQ{l,y\t) and FQp{y\t) 
satisfy the same functional equation, and therefore the coefficients of both series satisfy the 
same recursion. Since the first few terms of each coefficient sequence agree, this proves 
that these two classes, complete with the stopping height parameter y, are in bijection. That 
is 

Fs,(l;0 = //s(l,l;0. 

n 

Thus, the generating MGF HQ{l,y; t) and OGF //s(0 for half plane walks on a step 
set S are the same as the MGF F^^iy; t) and the OGF F^^it) for directed meanders on the 
projected set S^. There is one detail we have glossed over here: the x-coordinate of the 
steps. Because the half plane walks we are working with here are on sets of small steps, 
there is no problem with simply suppressing the x information. However, for more general 
sets of steps, the same holds. Since the only boundary is the .x;-axis, the y-coordinate of 
each step is what is important to allowing a walk to be 'legal' or not, and therefore the 
property that both classes are equinumerous will hold. Therefore, we may apply the results 
of Chapter |6] to show that "Kg has algebraic generating functions. 

7.3 Enumeration 

In view of our goal of enumeration, exact or asymptotic, and a corresponding generating 
function this is ideal: we obtain all the information we're after by considering the directed 
model found by setting x = I. The next example illustrates this for us. 

Example 7.2. Take as an example the step set shown in Figure ITTl S = {NW, NE, SW, SE}. 
A typical walk is shown in the same Figure. If we let H^ix, y; t) be the generating function 
for the half plane walks taken on S, then Hq will satisfy the functional equation 

HQ{x,y;t) = l+tixy + - + — + -]H,zix,y;t) - ti— + -|//s(x,0;0- 

\ y xy xj \xy y) 
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We take the horizontal projection by setting x = 1 to obtain 

Hs(l,y;t) = l+t(2y+-)HQ(l,y;t)-ti-\H^{hO;t). 

According to Proposition 17.2. ll this is the functional equation for the class of directed me- 
anders on f , the horizontal projection of S, allowing us to bring to bear the machinery 
given in Chapter |6l We suppress the x information and rearrange the equation to give 

H^(y; t) - t{2y + -)H.iy- 0=1- t-H^{0- 1), 

y y 

from which we may isolate the kernel in the denominator of the right hand side, like so: 

l-rf//s(0;0 

H^(y;t) = —^ —. 

1 - t(2y + j) 



We let 

K(y,t) = l-t(2y+-\ 

y 

and solve the characteristic equation K(y, t) = 0: 

y - f(2/ + 2) = ^ 2ty^-y + 2t = 0, 

1+ Vl - I6f 

^ ' = — s — • 

from which we take the small branch (the branch analytic at 0) with equation 

,, 1- Vl - 16?2 

y.it) = . 

We only need the enumeration information, so we ignore the parameter y by setting y = I. 
Then the generating function is 

"^^^'^''^-Y~^)- 4^(1^40 

from which we may find a formula for h(n) using the generalised binomial theorem. As in 
the directed case, this is a radical: 

4t-l+ Vl - 16^2 

«3W = J, . 
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Next, we need asymptotics. We may apply the black box Theorem l6.3.2[ First, we need 
the drift of the walk. It should be obvious that this is zero, but the general method is to 
find S = P'{1), so in our case 

P'(y) = 2 - 4 ^ ^'(1) = 

as expected. This tells us that we use the drift expansion for the asymptotics, given as 

4" 



[f ]Fe,(l; - vo^ f 1 + - + ^ + • • • ) 



where 



^ P"{X) 
This gives the asymptotic development for half plane walks on S as 

V2-4"/ c, C2 
-^ 1 + - + - 



[ni^e,(l;0 ^ 1 + - + ^ + 



The ability to exactly enumerate half plane walks on a step set S is a tool which we can 
use to find upper bounds on the exponential growth of the quarter plane walks on the same 
step set. After some background on the state of quarter plane walks. Chapter [8] will show 
how this is done. 



Chapter 8 

Walks confined to the quarter plane 



We have seen so far that walks in the whole plane have a rational generating function, 
easily found by an application of elementary analytic combinatorics. Walks in a half plane 
are solved using the bijection to a directed model which in tum is solved by the kernel 
method. These models have an algebraic generating function in every case. We finally tum 
our attention to the final region of the thesis: walks restricted to first quadrant of 7}. This 
class is much richer, and there are only the beginnings of a unified approach [il0l l6l. There 
are models with transcendental D-finite generating functions, and models with non-D-finite 
generating functions. 

In Section [8TT] we give the problem in analogy with the previous three cases, stating a 
functional equation that the MGF will satisfy. S ection W!2\ covers a reduction of the 1} - \ 
possible step sets to 79 non-equivalent models, classifying their associated generating func- 
tions as D-finite or non D-finite and providing an asymptotic growth template for the D- 
finite cases. As in the half plane model and the directed model, results will be dependent 
upon the drift. 

S ection [83] co vers the enumerative results for the D-finite cases, beginning with a survey 
of previous results. After the survey, we give in Section 18.3.21 proofs of the exponential 
growth for 8 D-finite step sets which form the base cases in our case analysis. Some of these 
proofs are original, and those taken from other authors are marked as such. Sections 18.3.31 
and 18.3.41 give an original systematic way of proving the exponential growth of each case 
by bounding the exponential growth above in a unified way and importing lower bounds 
from the base cases via a boot strapping method given in Lemmas 18.3.111 and 18.3.121 The 

55 
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results are summarised and notation for each case are introduced in Table 18. II 

8.1 Definitions 

Fix a step set S c ®. We call the class of walks taken on S restricted to the quarter 
plane Qq. It should be noted that on the step set S, we have Q^ c 'H^- This fact will prove 
useful in the asymptotic enumeration of the quarter plane walks. 

So, let q{n) be the number of walks on S confined to the quarter plane. Then the OGF 
for ^3 is 

n>0 

As with the half plane case, in order to find a functional equation, we need to parameterise 
the class. Let q(i, j; n) be the number of walks with steps from S, confined to the quar- 
ter plane. Then the MGF for Qq with parameters x and y marking x and y coordinates, 
respectively, of the endpoint, is 

The functional equation that Q<z{x,y\ t) satisfies is 



Qq{x, y;t) = \ + t{S {x, y)QQ{x, y; t) 

- {x^'']S{x,y)Q^{Q,y;t) 

- {y^''}S{x,y)Qs{x,Q-t)) 



+ 6:>(S)— es(0,0;0, (8.1) 

xy 



where 

a>(S) 



1 if (-1,-1)6 6, 

otherwise. 

Note the diff'erences from the half plane functional equation. The first diff'erence to note 
is the term t[x'^''^}S {x,y)Q(z{0,y;t), removing the walks of length n + \ which make their 
way into the region x < 0. We also have the term a)(S)— 2s(0, 0; f), which corrects for 
undercounting if the SW step is included in S and therefore removed twice since it is 
both x-negative and j-negative. 
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Again, we have a notion of drift of a step set, this time with two components. We recall 
it in the following definition. 

Definition 19. For a given step set S, the drift is defined to be the vector sum of the steps 

in S. 

Note that since the drift is a sum of 2-vectors, it will be a 2-vector also. Thus we give 
the drift in terms of the sign of each component. This is illustrated by an example. 

Example 8.1. Take the step set n". The horizontal and vertical steps will have vector 
sum 0, so the vector sum of the whole step set is in the NW direction, given by the ordered 
pair (-1,1) or the vector diagram n>. In this case, we say the drift is j-positive and x- 
negative. 

This new functional equation is a little more complicated than before. First of all, it's 
in terms of at least three generating functions and in some cases four. The success of the 
kernel method in the directed and half plane models is due to the low number of unknowns 
that can be made dependent by solving the kemel equation, a property missing from this 
case. We can rearrange the functional equation to isolate QQ(x,y; t), giving 

1 - t({x<''}S(x,y)QdO,y;t) + {y<'}S(x,y)Qs(x,0;t)) + a;(S)^2e(0,0;0 

Qs(x,y;t) = — — ■ . 

l-tS(x,y) 

(8.2) 
In analogy with the directed model, the denominator is called the kemel of the model, 
and plays a part in the orbit-sum enumerative techniques used by Bousquet-Melou and 
Mishna in ^. The following section gives an overview of the part the kemel plays in the 
classification of Qzix, y; t) as algebraic, D-finite or non D-finite. 

8.2 Classification 

In the full plane model, there are 2^ - 1 possible step sets S available. Of these, some 
will give rise to trivial classes Qq. These are subsets of {SE, S, SW, W, NW}. There are 
also step sets which will never interact with one boundary, making them automatically half 
plane models. Shown in Figure [STI is one of the largest step sets that will behave in this way. 
Another is found by reflecting the first across j = x. By removing trivial step sets as well 
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as step sets which are automatically half plane models, thus solvable by the kemel method, 
and considering symmetries, Bousquet-Melou and Mishna in [lZjI reduce this number to 79 
non-equivalent true quarter plane models. 



K 



Figure 8.1: The largest step set which is automatically a half plane model. 

Associated with each model is G(S), some dihedral group of birational transformations 
which preserves the kemel K{x, y;t) = l-tS (x, y) of a model - hence dependent on S. Note 
that not all dihedral symmetries are in G(S), the group depends on the generators defined 
in [|7l. For 23 step sets S, this group is finite while for the remaining 56 it is infinite. 
There is a lot of evidence for a correspondence between the cardinality of G(S) and the 
classification of the associated generating function: it is shown in [|71 that 22 of the 23 finite 
group models have D-finite generating functions - with the 23rd requiring another approach 
in its own paper [|5l - and it is shown in [6j that 51 of the 56 infinite group models have non 
D-finite generating functions. The remaining 5 models are known as the singular models, 
and are particularly stubborn when it comes to being shown they are non D-finite (see [|22]| '). 

For the remainder of this chapter, we focus on those models associated with a finite 
group or those with a D-finite generating function. These models are listed in Table 18.21 
along with experimentally found asymptotic expressions, first given in flU, which we shall 
discuss in the next section. For the remainder of this section, we give a stronger classifica- 
tion of the generating functions, which helps with the asymptotic enumeration later on. We 
recall a definition from dU. 

Definition 20. flU A power series Q(t) = 2n>o '^nt" in QE/^1 is called a G-series if: 

1 . it is D-finite; 

2. its radius of convergence in C[[f]] is positive; 

3. there exists a constant C > such that for all n 6 N the common denominator of 
ao, ..., fl„ is bounded by C". 
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Since generating functions are elements of Z[[f], they trivially pass the third condition. 
Furthermore, the generating functions in question here are a priori D-finite, so the first 
condition is satisfied also. We can show the second condition to be satisfied by considering 
the final line of the previous paragraph: the number of walks of length n cannot exceed | S|". 
We summarise this in the following proposition. 

Proposition 8.2.1. /@ Theorem 1] For a given step set ®, the radii of convergence of each 
ofQe(0,0;t), Qe(0, Ut), 2s(l,0;0, 2s(l, 1;0 is bounded below by j^. 

Proof. The total number of walks in "W taken on S is |S|", thus this forms an upper bound 
on the number of walks in Qq taken on S, which is also the coefficient q{n). By virtue of 
the exponential growth formula from Theorem 14.3.21 this upper bound on the coefficients 
of the generating function forms a lower bound on the radius of convergence, proving the 
statement. n 

This satisfies the second property of G-series, completing the proof of the following 
theorem from [|4j. 

Theorem 8.2.2. /@ Theorem 2] The generating functions for the quarter plane models Qq. 
for i = 1, ..., 23 are G-series. 

In im, Bostan and Kauers show Theorem l8.2.2l along with a property of the asymptotic 
growth of the coefficients of the generating functions. It can be concluded from Theorem 
18.2.21 and results in [f T2] that the coefficients of any one of the generating functions in 
question fits the asymptotic template 

^(n)~/<-y8V(lognf. 

Furthermore, Bostan and Kauers give compelling evidence that y = for each of the cases 
we consider in this chapter, thus these models probably have the same growth template as 
an algebraic series. We discuss this further in the following section on enumeration. 

8.3 Enumeration 

We begin the enumeration discussion by surveying some of the more general previous enu- 
merative results linked with the classification of the generating functions in Section [8.3. II 
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These are split into rigourous results given by Bousquet-Melou and Mishna in [|7]| and the 
experimental results of Bostan and Kauers in [j4l. 

Sections [8.3.2l through 18.3.41 contain the original contribution of this thesis, where we 
prove the exponential growth factors for the D-finite models by bounding the exponential 
growth and applying a squeezing argument. Upper bounds are found in a unified way via a 
relaxation of the restrictions placed on the class. Lower bounds are given by a case analysis 
reduced from 23 to 8 cases thanks to Lemmas [8.3. Ill and l8. 3. 121 

8.3.1 Previous Results 

In this section, we give some rigourous and experimental results of previous authors, high- 
lighting how they are tied to the classification of the generating functions. 

Rigourous results using orbit sums 

In Chapter|6l we used the kernel method, which takes a solution to the kemel equation and 
uses it to bind t and the catalytic variable y, which marks the altitude of the endpoint, in 
such a way that we reduced the number of unknowns and found a solution. For 22 of the 23 
finite group models, Bousquet-Melou and Mishna apply a similar process to the functional 
equation satisfied by the generating functions of those models in [|71. We briefly discuss the 
process here, leaving the details to the source material. 
First, we rewrite Equation (18.21 ) as 

K(x,y;t)xyQz(x,y;t) = xy - xyt[{x''^}S(x,y)Qs(0,y;t) 

+ {y'''}Six,y)Qs(x,0;t)) + tco(Q)Q5iO,0;t). 

The multiplication by xy will clear x and y from the denominators in every term, thus 
xy{x'^^}S(x,y) is now a polynomial in y and similarly xy{y^^}S{x,y) is a polynomial in x. 
Now let F{x) = xy{y^^}S(x,y)Qz(x,0; t) and G{y) = ;cj{x<°}5(x, j)2e(0, j; t), so the above 
equation becomes 

K{x, y; t)xyQ^ix, y; t) = xy - tF(x) - tGiy) + co{<Z)tQs{0, 0; t). (8.3) 

Because the unknowns are no longer mixed as functions of x and y, the resulting equa- 
tion is well suited to the orbit sum method. 
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The orbit sum method 

The group G(S) is defined as a group of birational transformations of the plane preserving 
the kernel K{x,y; t) = \ - tS{x,y). This is a dihedral group with two generators 

O : {x,y) ^ (x,y) and ¥ : {x,y) ^ (x,y'), 

where the terms x' and y' are rational functions in x and y and depend on the step set S. 
The fact that each generator acts only on one of the variables x and y is the key concept 



here. Acting on Equation (|8.3I) with O, we find 

K(x, y; t)x'yQs(x',y; t) = x'y - tF(x') - tG(y) + a)(Q)tQs(0, 0; t). (8.4) 

Then we take the difference of Equations (18.31) and (18.41 ). eliminating the term tG(y): 

K{x, y\ t) {xyQ^ix, y\ t) - x'yQ^ix ,y\ t)) = xy - x'y - tF{x) + tF{x'). (8.5) 

The next step is to act on Equation (18.41) with *F, generating a new equation, and adding it to 
Equation (18.51) . continuing to alternately act with O and *F and taking the altemating sum. 
If the group is finite, as it is for our 23 cases, this process will terminate and all unknown 
functions on the right hand side will vanish. This leads to the equation 

y sgn(g)g • {xyQs{x,y; t)) = — V sgn(g)^ • {xy). 

The coefficient of each summand is sgn(g) = (-1)^, where k is the number of generators O 
and *P in the decomposition of g. For 19 cases, this is enough to show D-finiteness. Es- 
sentially, the generating function is expressed as the constant term of a multivariate rational 
function. Since D-finite functions are closed under Hadamard products 11201 . it follows 
quickly. For enumerative purposes, we can arrive at an expression for the generating func- 
tion as a rational series extraction 



xyQs{x,y\ t) = {x^^} ^ V sgn(g)g ■ (^3;). 

Kix,y;t)^ 

From this, we may derive the formula for the number of walks of length n ending at (i, j): 



q{i,j;n) = [x'^y^'] 



i+i,,;+i- 



Yj sgn(g)g • (xy) 



\gsG 



S{x,yr. 
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For the remaining four models the approach isn't as useful because the right hand side 
for each vanishes, leaving 

geG 

For one of these models known as Gessel's model, 623^ we have not mastered the use 
of rotational symmetry and so this model cannot be handled by this method in a straight 
forward manner. Bostan and Kauers treat this model in [|3. For the remaining three, it is 
due to an x/y symmetry: the presence of {y, x) in the orbit of ix,y). Bousquet-Melou and 
Mishna approach this with an alternating sum over half of the orbit instead, and manage to 
solve these three models. 

While these formulae and expressions provide invaluable tools to generate enumerative 
data and a classification of the generating functions for 22 of the 23 cases, the aim was not 
to provide asymptotic data. The rational series extraction gives no hint as to the location of 
the singularities of the generating function. Series generation is easy to do with help from 
computers and this leads to another approach. The next section will outline a very succesful 
experimental regime employed by Bostan and Kauers to this end. 

Bostan and Kauers' experimental approach 

Recall that the asymptotic growth of a counting sequence is in correspondence with the 
radius of convergence of the associated generating function. In [4J, Bostan and Kauers 
use the D-finiteness, as well as other readily available information about the model, to 
conjecture the radius of convergence of each model and the asymptotic expressions. 

Recall that a D-finite generating function satisfies a linear differential equation of or- 
der r with polynomial coefficients c,(0, where Co(0 is the leading coefficient. If we assume 
the existence of a simply connected domain on which the coefficients are analytic, the ex- 
istence theorem on page 3 of [|26l guarantees that in a neighbourhood of a point to such 
that co(?o) i^ 0, there exist at most r linearly independent analytic solutions to the above 
differential equation. This implies that singularities can only occur at points which are 
roots of the leading coefficient Co(t), giving an upper bound of deg(co(0) candidates for the 
dominant singularity. 

Next, recall that Proposition l8.2.2l bounds the radius of convergence of a G-series below 
by 1 /I S|, further reducing the number of candidates for the dominant singularity. 
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So, to use these facts Bostan and Kauers required a differential equation satisfied by the 
generating function of a given model. By generating all walks on a given step set up to 
length N (at that point, the maximum achievable was around N = 1000, but now with the 
help of Steve Melczer, it is around A^ = 2000) and searching the Weyl algebra Q[t]{Dt) of 
differential operators satisfied by the series modulo t^, they were able to find a differential 
equation satisfied by the partial sum. The idea is that if the generating function happens 
to be D-finite, then for a sufficiently large A^, a differential equation satisfied modulo t'^ 
will provide a differential equation which is really satisfied by the full power series. Anal- 
ysis of the polynomial coefficients of the found differential equation yields the asymptotic 
information found in Table 18.21 

The column of exponential growth factors /3 provide us with strong candidates for the 
rigorous results. In fact, while these results are not proved rigorously in ^, the approach 
to do so is similar to that used in a paper by the same authors on Gessel's walks [ISl. They 
are omitted from [l4l because the computations required are quite large, motivating the 
development of another approach with more elementary techniques. 

Our approach, developed over the next three sections, uses restrictions and relaxations 
of the constraints on a class to bound the number of walks by known quantities, and a small 
amount of analysis to squeeze the exponential growth. Thanks to Lemmas [8.3.11l and l8.3.12l 
in Section [8". 3. 41 we are able to reduce the cases analysis required to 8 base cases. We begin 
by proving the exponential growth of these base cases in the following section. 

8.3.2 Base cases 

In Section 18.3.41 we give a way of importing lower bounds on the exponential growth of a 
step set S from proper subsets of S. This gives us a way of reducing our 23 cases down to 8 
base cases which are true quarter plane walks. There are two models which each import 
exponential growth from a half plane model, but in both cases the model is a class of Dyck 
prefixes, so we don't prove these base cases here. This section will prove the exponential 
growth factors for each quarter plane base case, first giving the base cases we can fully 
enumerate, then giving those for which we must enumerate a sufficiently large subset. 
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Table 8.2: The number of walks ^,(«) is conjectured to grow like K/S'^n" ([4, Table 1]). 
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Full counts 

There are two base case step sets which have quarter plane models susceptible to direct 
counting techniques. These are S2 and S^, pictured below in Figure |8^ 



Figure 8.2: The step sets which we count directly. 

Let us move through them in order. The first, S2 = X is pretty straightforward. 

Proposition 8.3.1. The exponential growth factor for the counting sequence q2{n) is 4. 

Proof. For S2, we have a simpler count. If we take a walk on S2 and look at the projection 
onto each axis, we see that the endpoint of the walk traces out a Dyck prefix on either axis. 
Moreover, since every step is positive or negative in both x and y, the endpoint projection 
moves on both axes at every step. So, if q2{n) is the number of walks of length n on S2 
and d{n) is the number of Dyck prefixes, then 
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gives the result. 









Finally, we move on to Sn = A. The proof is taken from an informal discussion in [|7], 
which can also be found in [[211 Section 3]. 

Proposition 8.3.2. ^0 Section 5.2] The exponential growth factor for the counting sequence 
qn(n) of walks on S^ = A is 3. 

Proof. To count walks on Sn = A, note that before we may use W, which is x-negative, 
we must use an .x;-positive step. But the only x-positive step is SE, which is y-negative, and 
before we may use it we must use the N step. Thus, the number of N steps must outweigh 
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the number of SE steps which must in turn outweigh the number of W steps. This puts 
the walks in bijection with Standard Young tableaux of height at most 3, which we explain 
presently. 

Standard Young tableaux are labelled Ferrer's diagrams of a partition such that the boxes 
are labelled in a strictly increasing manner from left to right and from bottom to top. We 
construct a standard Young tableau of height three of size n from a walk w = (si, S2, ..., s„) 
as follows. If Si = N (resp. SE, W), then place label i in the next available space to the right 
on the bottom row (resp. second, top). The final tableau is increasing from left to right by 
construction, and the prefix condition #N > #SE > #W ensures that it is increasing along 
the columns. 

Standard Young tableaux are well understood and may be enumerated via the hook- 
length formula. The number of standard Young tableaux of height 3 and size n is the nth 
Motzkin number. An application of Stirling's approximation will give the result. n 

Next, we move on to the longer list of cases for which we must count a sufficiently large 
subset. 

Counting a subset 

For the majority of the base case step sets, closed form enumeration of all walks of a given 
length is beyond the reach of current techniques. We look instead for easier to count subsets 
with the same exponential growth. Some will be walks returning to the origin, as in step set 
1 1 , others will be walks ending on a boundary. Often, we appeal to an upper bound given by 
the half plane walks taken on the same step set. For the proof that these upper bounds have 
the same exponential growth, see Theorem 18. 3. lOl in Section 18.3.31 We begin by counting a 
subset of S5 = {NE, S, NW} by injecting walks on another step set. This isn't really a base 
case, but the method of proof is difi"erent enough that it warrants inclusion here. 

Proposition 8.3.3. The counting sequence qsin) of walks on S5 = \ has exponential 
growth factor 3. 

Proof. Note that we may not count walks returning to the origin: these have the same count 
as walks taken on the negative drift counterpart of this steps set, S15, so we must be a little 
clever. Consider for a moment the step set T = {NW, E, S}, shown in Figure [83l This set 



CHAPTERS. WALKS CONFINED TO THE QUARTER PLANE 68 



Figure 8.3: The helper set 

is a reflection of Sn across the line y = x, so the number of walks of length n is the nth 
Motzkin number. 

Note that if we rotate a step so it takes us further from any boundary we are doing no 
harm. So if we rotate the E step so it becomes a NE step, all allowed walks from 7~ are still 
allowed on the new step set. This is the same as injecting the walks taken on 7~ as a subset 
of the walks taken on S5 with the function 

f:T ^ Qs 

NW ^ NW 

E ^ NE 

S ^ S. 

Thus, we have 

Mn < qsin) < hsin). 

Taking lim ^ log of each term gives 

n 

Moving on to the next base case, Sn = {N, SE, S, SW}, which is the first negative drift 
base case. 

Proposition 8.3.4. The exponential growth factor of the counting sequence qi\{n)for quar- 
ter plane walks on Q\i = Mis! y3. 

Proof Consider the walks on this step set which return to the origin. Notice that any walk 
returning to the x-axis must be of even length, since every step is either y-positive or y- 
negative. Thus, any walk returning to the origin must also be of even length. We begin by 
considering a walk of length 2n taken only on the N and S steps, returning to the origin. 
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This is simply a projection of some Dyck path onto the y-axis, and these are counted by the 
Catalan numbers. 

We replace some of the S steps by SW and SE steps. Since the walk must return to the 
origin, the walk on {SW, SE} must form a Dyck path with respect to the y-axis, which must 
also be of even length. These too are counted by the Catalan numbers. Moreover, we must 
choose some even number of the j-negative steps to replace, so there will be a summation 
of binomials also. Explicitly, the number of walks returning to the origin on Si i is given by 

[n/2i , ^ . 

^n(0,0;2R) = C„_^rb. 
The quantity 

Ln/2J / ,, 

is an expression for the nth Motizkin number M„, so we may rewrite the above as 

^ii(0,0;2n) = C„M„. 
By applying Stirling's approximation, we know that for C, D,r,s£ R, we have 

M„ ~ C3V and C„ ~ D4V, 
so given some < 6 < 1, we know there exists some A^ sufficiently large such that 

CD(2 V3)2"n'-+-^ 

for all n > N, that is, the exponential growth factor of ^ii(0, 0; 2n) is 2 VS. 
Thus 

^ii(0,0;2n)x(2V3)^ 

and we may apply lim ^ log to each term in the inequality 

^ii(0,0;n) < qu(n) < hu{n), 
taking care to take the limit supremum on the left, to prove the result. n 

For Si3 = X, another negative drift step set, we look to the projection of the walk on to 
either axis, similar to the argument for 62^ but with a bit more analysis. 
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Figure 8.4: A walk taken on S13 shown with the projection onto each axis. 



Proposition 8.3.5. The counting sequence qij(n)for QP walks on S13 = X has exponential 
growth factor 2 V6. 

Proof. As in S2, we look at the projection of a typical walk on the x and y axes. Since all 
steps are either 3;-positive or negative, the projection onto the y axis forms a Dyck prefix. Let 
us consider the projection onto the x-axis. There is one step with no horizontal component, 
and so the projection will form a Motzkin prefix, but with the stationary steps forced to cor- 
repsond to the downward steps in the projected Dyck prefix. Figure [84l shows an example 
of a walk returning to the x-axis with the steps in blue highlighting this correpsondence. 

In order to simplify computations, we consider the walks ending on the .x;-axis. Thus 
the projection onto the y-axis forms a Dyck path, something that we are quite familiar with 
counting: let the number of projections onto the j-axis be Py{n). Then Py{n) = C„. The pro- 
jection onto the x-axis forms a Motzkin prefix, and as already mentioned, the 'stationary' 
steps correspond only to y-negative steps in the original walk. 

Now count the number of Motzkin prefixes with this restriction. The walks return to the 
X-axis, and since all steps are either y-positive or negative the walks must be of even length. 
Of the n downward steps, k are chosen to be the steps with no x component. The remainder 
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of the 2n steps must form a Dyck prefix, giving the formula 

Px(n) = ^i \di2n-k). 

From Chapter |6l we know the asymptotic growth for Dyck prefixes is 

d(2n -k)-- C2^"-\2n - kf, 

where C is a positive real constant and r < 0. Thus for < 6 < 1, there is some A^ such 
that n > N implies 

N 



J] hd{2n -k) + j] ( j(l - e)C2'"~'{2n - kf 

k=0 ^ ' k=N+l ^ ' 

< pM < £ (!)^(2n -k)+ J] h(l + e)C2'"-\2n - kf. 

k=0 ^ ' k=N+\ ^ ' 

Let's look at the lower bound first. We have 
pM > J]hd(2n-k)+ J]h(l-e)C2'"-'(2n-ky, 

k=Q ^ ' k=N+\ ^ ' 

> Yj (ik^n -k) + {l- e)C2"\2ny J^ (?) i\] since r < 0, 

k=0 ^ ' k=N+l W \^/ 

- Zi ft) r*-^" ~^''~ <^22"(2n)'" j -) + (1 - e)C2^"(2nYi-\ binomial theorem, 

- Z (^) (^^2'^ - ^^ - C2^"(2ny (^) ) + (1 - ^)C( ^f"i2ny. 

Now, the sum with upper bound of summation N is the error term. Note that the constant 
upper limit of summation makes the growth as n ^ oo of the positive portion polynomial, 
and the negative portion has its exponential growth dominated by ( V6)^". Thus there exists 
some new A^~ such that n > N~ forces 

Y, (^) [d{2n -k)- C2^\2ny |i j j < |c( <6f\2ny , 
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so for n> N~ 

{\-el2)C{<6f\2nr<pM)- 
The upper bound is similar, giving 

pAn)<{l-e)C{<6f" 

for n sufficiently large. This proves 

pM)x{<ef". 

Recall that the number of walks returning to the ;c-axis on S13 is the product 

2_^ qi3(i, 0; 2n) = py(n)px(n). 

i>0 

We already know that Pyin) ~ ^3^^" , and we just showed that Px(n) ~ D( V6)^"(2n)' for 
some positive D and s <0,so for some D' and s' we have 



^^i30',0;2n)~D'(2V6f«(2n/. 



i>0 



Thus given the inequality 

2_j ^13(0, 1; n) < qi3(n) < hu{n), 

!>0 

we know that for n sufficiently large we substitute each quantities respective growth func- 
tion, and working in the usual way (taking the lim sup on the left) we find 

ySi3=2V6. 

n 

This brings us to S15 = A. We find the exponential growth by considering walks 
returning to the origin. 

Proposition 8.3.6. The exponential growth factor for the counting sequence qi5(n)for QP 
walks on Q = J\ is 2 \2. 
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Proof. The walks returning to the origin are easy to count. Every walk which returns to the 
origin must be of length 4n, since there must be the same number of y-positive steps as y- 
negative, and the same number of jc-positive as ;c-negative. Now, if we look at the projection 
of a walk that returns to the origin onto the j-axis, the endpoint traces out a Dyck path of 
length 4n. Then we look at the subwalk taken on the x-negative steps. This is a Dyck path 
of length 2n. So, we have 

^15(0, 0;4n) = CinCn - C4^"r(4ny = C(2 ^f2f"(4ny 

for some C 6 R^ and r < 0. We then get 

^15(0, 0;n) < qisin) < hi^in), 

which we may work on in the same way as number Sn: by taking the limsup on the left, 
and getting 

ySi5=2V2. 

n 

The next base case is S19 = {NE, S, W}. This is a very famous model of a quarter plane 
walk, and is known as Kreweras' model, since he was the first to study it. We prove this 
factor along with the factor for S20, since one step set is the reverse of the other. 

Proposition 8.3.7. / fJTl/ The exponential growth factor of the counting sequences qi9(n) 
and q2o{n), for QP walks on S19 = \ and S20 = /*" respectively, is /3i<) = /32o = 3. 

Proof. In 1965, Kreweras proved in [17J that the number of walks on S19 of length 3n, 
terminating at the origin, is given by 

- 4" (3n\ 

" ~ (n + l)i2n+ 1)\ n j' 

An application of Stirling's approximation tells us that the exponential growth of this quan- 
tity is 3^", as desired. This provides a lower bound on the exponential growth of the unre- 
stricted quarter plane model, giving the first part of the result. 

The next base case to count is the step set S20 = {N, E, SW}. Note that this is the 
previous case in reverse, so the number of walks returning to the origin of length 3n is k„, 
which we just showed has an exponential growth factor fijo = 3, giving a lower bound on 
the quarter plane model and proving the second part of the result. n 
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The next base case, S22 = {E, SE, W, NW} is an interesting case. It is known as the 
Gouyou-Beauchamps model. There are a few different ways to obtain the lower bound we 
need here. The proof is taken from [jTJ. 

Proposition 8.3.8. /0 Proposition 11] The counting sequence qiiin) ofQP walks on S22 = 
"^ has exponential growth factor ^22 = 4. 

Proof. First of all, we notice a bijection between this model and pairs of non-intersecting 
Dyck prefixes. To pass from a pair to a quadrant walk, we must assign a direction to each 
possible pair of steps. Hence, the number of walks ending at a given position {i, j) is given 
as a 2-by-2 Gessel-Viennot determinant, and we can expect a closed form expression for 
this number. 

Next, note that we may make the transformation (/, j) \-^ (i + j, j) to take quarter plane 
walks injectively to walks in Si = {N, S, E, W} lying between the jc-axis and the line y = x. 
This is the form they were studied in by Gouyou-Beauchamps, who provided a formula 
for the number of n step walks ending on the x-axis as a product of Catalan numbers 
in [|T3]| . This interest grew from a far from obvious bijection between these walks and 
Young tableuax of height at most 4. The Gouyou-Beauchamps closed form expression for 
walks ending at the origin is 

^(0 0-2m) = 6(2m)!(2m + 2)! ^ 24-4^'" 

m\{m + iy.{m + 2)l{m + 3)1 nnv" 

This has the desired exponential growth, and so may be used as a lower bound to prove ;022 = 
4. n 

This concludes the proof of the exponential growth of the base cases. In the following 
section, we show how to tightly bound above the exponential growth of each step set by 
relaxing the regional constraints on the class. 

8.3.3 Upper bounds on the exponential growth 

We wish to use the lim - log technique to squeeze the exponential growth of the counting 
sequences for each S,. In this section, the main result is the following. 

Theorem 8.3.9. (Upper bound) For i = 1,...,23, the exponential growth factor for the 
number of walks taken on ®, is bounded above by the value Pi given in Table WT\ 
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Figure 8.5: Removing the y-axis as a boundary to produce an upper bound on the exponen- 
tial growth 



Now, we know intuitively (and we proved in Section lOI) that since the number of walks 
on a step set S cannot exceed |S|", if jS is the exponential growth factor for the number of 
walks on S, we have 

;0<|S|. 

It turns out that for most of the 23 D-finite cases this bound is tight. We were intrigued 
by the fact that for those S, with 3^-negative drift, this bound is too large. This is because 
bounding by the number of steps is the same as removing both boundaries, and considering 
(3s as a subset of 'W^, and the enumeration of 'Wz doesn't depend on the drift, and as 
can be seen from Chapters |6] and Ul when a boundary is introduced drift influences the 
enumerative properties of a model. 

As an alternative bound, we consider removing only one boundary, as shown in Figure 
18. 5[ Note that for i e {1, ..., 23}, if the drift is non-zero it is in the j-direction only and 
therefore will be preserved by a horizontal projection. Thus, considering Qq as a subset of 
'Hq will give drift dependent upper bounds on the exponential growth. 

For a given step set S,, viewing a quarter plane walk as a restricted half plane walk will 
give the inequality 

q(n) < h(n), 

where q(n) and h(n) count quarter and half plane walks, respectively, of length n. Then, for 
the right C, D,/3, 7, r, s, we have 



q(n) ~ Cfi"n' and h(n) ~ Dy"^^ 
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and for sufficiently large n we will have 

C/3V < DyV. 

Once again, taking lim ^ log on both sides will yield 

logyS < log y ^ P <y. 

Now recall that the exponential growth factor for half -plane walks with zero or positive drift 
is equal to the cardinality of the step set, so this bound is as good as the previous one for 
these cases. Also, we expect it to improve in the negative drift cases since for the class 'Hq 
the drift afi'ects the asymptotic growth. Let us consider an example. 

Example 8.2. We showed in Proposition 18.3.41 that Sn = /k has a lower bound of 2 VS 
on the exponential growth, and appealed to an upper bound given by the half plane walks 
on the same step set. Here we produce this upper bound by appealing to the theorems on 
directed walks from Chapter |6l applied through the bijection given in Chapter |7] to half 
plane walks. 

First, we represent Sn by its inventory 

X 1 1 

S{x,y) =3; + - + - + —, 

y y xy 

and take the horizontal projection by setting jc = 1, taking the directed inventory. This gives 
us the weighted inventory for the directed meander on the horizontal projection V, 

P(y)=y + -. 

y 

We find the drift to be negative: 

P'iy) = 1 - 4 ^ P'iD = -2 < 0; 

y^ 

so Theorem 16.3.21 tells us that the exponential growth factor of hn(n) is equal to P(t), 
where r is the positive solution to the equation 

P'(T) = 0. 

So, we let 

1-4=0 

y^ 
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which is easily solved for 

T= Vs. 

Plugging this back into P(y) gives 

y = P(T) = 2^f3, 
which is tight according to Proposition 18 .3 .41 

So in Example 18.21 using the relaxation to a half plane model, we have found an upper 
bound for jSii which is an improvement on the upper bound provided by the number of 
steps. It is an elementary application of the method given in Example 1 8. 2 1 to the remaining 5 
negative drift cases to verify that it is a tight upper bound in all cases. This leads to the 
following result, a restatement of Theorem 18 . 3 . 9 1 from the beginning of Section [8.3.3l using 
a little more precision. 

Theorem 8.3.10. The exponential growth factor of a quarter plane walk with a set of small 
steps S,-, / = 1, 2, 3, ..., 23, is bounded above by 

L the number of steps for a walk with non-negative drift, 

2. the quantity P{T)for a walk with negative drift, 

where P(y) = S(l,y)is the inventory of the horizontal projection of S,- and r is the structural 
constant of the half plane model. 

The proof of Theorem l8.3.10l is simply an application of the method of Example 18.21 in 
each of the cases and so is omitted. Now that we have a method for providing upper bounds 
which agree with the findings of Bostan and Kauers in all cases, we move on to finding 
lower bounds. This is a little more complicated than the upper bounds, and is detailed in 
the following section. 

8.3.4 Lower bounds on the exponential growth 

In this section, we show how to find lower bounds on the exponential growth for the given 
models. In some cases, we are able to use a restriction, such as the walks returning to the 
origin or those returning to an axis, but this doesn't work for all cases. We are then forced 
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to construct lower bounds from step subsets that we know the exponential growth of. We 
illustrate the intuition of this in the following example. 

Example 8.3. Take Q5 = Y, which we shall call the positive drift fork. Notice that if we 
take any walk of a given length on S5 and insert a north step at any internal point, the 
remainder of the walk moves no closer to the y-axis and further from the x axis, giving us a 
walk from S7 = 1. If ^5(2) is the number of walks of length / on S5, and ^5(/) ~ C3'i\ we 
can build walks of length n in S7 by taking walks of length i on S5 and inserting n - i north 
steps. See Figure [8l6] for an example, where we insert N steps at the positions marked by a 
blue circle. The number of walks on S7 built in such a way is 






where we denote it by q' since in general, this method might not give all possible walks on 
the larger step set. Now, intuitively, since q^in) has dominant growth 3", we would like to 
get q'lin) with dominant growth 4" by applying the binomial theorem, like so 
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Figure 8.6: Inserting north steps at points of a walk on S5 to create a walk on S7. 



Furthermore, we could take S7, the trident, and increase the number of steps by two, 
inserting a pair of steps in opposite directions. The only thing we need to make it work 
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is that the subwalk taken on the new steps is a meander with respect to the relevant axes. 
If d(i) is the number of Dyck prefixes of length i, and we are adding the steps {E,W} to 
the step set S?, we get q'^(n), the number of legal walks from Sg = {N, NE, E, S, W, NW} 
built in this way, counted by 



,->n V/ 



r>0 

Since qjin) has exponential growth 4" and d{n) has exponential growth 2", we'd like to 
replace the functions with their dominant growth and apply the binomial theorem, as before, 
to get a set with dominant growth 6". Again, doing this requires some care. Figure 18.71 
shows an example of this process, where we insert horizontal steps at the positions marked 
by a blue circle in the graphic on the left. Note that the subwalk taken on the blue steps is a 
meander. 





Figure 8.7: Inserting a meander on {E, W} into a walk from £7 to obtain a walk on Sg. 



The intuitive idea has merit, but the asymptotic expressions contain a subexponential 
factor. We show the details in Lemmas 18.3.111 and 18.3.121 The following section outlines 
how restrictions on the stopping conditions for a path work for some step sets, but can fail 
for others. 
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Restricting stopping conditions 

Given that upper bounds were obtained in a unified way by considering the relaxation of the 
problem to a half plane model, a good first idea for finding lower bounds would be finding a 
restriction with the same exponential growth. For example, this was the technique we used 
in Proposition 18.3.41 to prove the exponential growth for Sn. Unfortunately, restricting the 
stopping condition fails to provide tight lower bounds for the positive drift cases S,, i e 
{5, 6, 7, 8, 9, 10}. This is because each positive drift case is the reflection of a corresponding 
negative drift case, as shown in Figure [ 



^ — ^ 



Figure 8.8: Sy and S 



11- 



What this symmetry implies is a bijection between the walks returning to the origin on 
either step set. Following a walk ending at the origin on one step set in reverse will give 
you a walk returning to the origin on the other. Given this bijection, the exponential growth 
factors must be the same, and we cannot use walks returning to the origin for a tight lower 
bound on the exponential growth factor of step set 7. 

In fact, this is true of every pair of step sets which are reflections across a horizontal 
axis of symmetry. To see this, we note that walks returning to the origin are in bijection 
for every pair, and so the exponential growth factor for such walks is bounded above by the 
exponential growth factor for general quarter plane walks on the negative drift version of 
each pair. We already showed in section [03] that these factors are in turn bounded above 
by the exponential growth factors of the number of half plane walks on the same step set, 
each of which is strictly smaller than the number of steps. Thus, for every positive drift step 
set, we need a better lower bound than the number of walks retuming to the origin. The 
following section will provide these. 

Lower bounds from subsets 

Recall Example 1 8. 3 [ By carefully inserting new steps (see Figures [831 and [877]) into walks 
on a step set, we were able to bootstrap a lower bound on the exponential growth of a larger 
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step set. This relied on the models having only exponential growth, so the proof reduced 
to an application of the binomial theorem. In our 23 examples, the growth has a constant 
factor and a subexponential factor. While the analysis does get a little messier the relevant 
part of the result is still the same. 

Lemma 8.3.11. Let q(i) ~ CjS'f, where s <0,C e R'*',/3 e [1, oo), and define 



i>0 ^^' 



Then q'{n) has exponential growth factor fi + 1. 



Proof. Fix 6 > 0. Then since q{n) ~ Cfi"n'\ we know that there exists some A^ such 
that n > N gives 

We take n > N, and consider the quantity 



C/3"n' 



q'(n) = J^hq(i). 

We may decompose the summation into two parts, those terms with summation index 
smaller than A^^, and those larger. We look first at the lower bound given by the growth 
properties: 

q'(n) > J]hq{i)+ J]i%l-e)C0n, 

> Yj (")^(') + (1 - ^)<^"' Yj ("r ''^^^^^ "' - ''' 

i=0 ^^' i=N+l ^^' 

(=0 ^^' i=0 ^^' 

N I V 

- Yj r K^(') ~ ^"'^') + (1 ~ ^)^^ + ^^""' (^y *^ binomial theorem). 

We must now show that the error terms become insignificant as n increases. We do 
so similarly to before: we divide throughout by C(J3 + l)"n^ and take the limit. Since the 
upper limit of summation is a constant, the growth of the numerator is polynomial, and thus 
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dominated by the growth of (yS + 1)". That is, we may find anN > N such that forn>N 
we have 



j:lo("i)(q(i)-Cn'0) 



Cifi + lyn' 



e 
< -. 

~ 2 



So taking n > N , we find 



(l-26)<(l--)< 



q'(n) 



We work similarly on the upper bound. The main difference is that we wish to increase 
the value of the bound when making simplifications, so rather than replacing the polynomial 
factors with something smaller, we remove them entirely. 

q'{n) < £["W)+ ^ h((l+6)C;e'n, 

M) ^^' i=N+\ ^^' 

(=0 ^ ' i=N+l ^ ' 

N I s 

1=0 ^^' 

We then use the same analysis as before to show that the error term becomes smaller than 
|C(jS + 1)". That is, there is some A^"^ such that taking n > N^ gives 



q'{n) 



<(l + y)<(l+26). 



ifi + 1)" 

Putting it all together, we have that for n > max {A^^, A^ }, the quantity q'{n) satisfies the 
inequality 

(1 - 2e)C(fi + \fn' < q'{n) < (1 + 2e)C(fi +1)". 

Consider what happens when we take the logarithm of the bounds, divide by n and take the 
limit as n tends to infinity, our standard method of isolating the exponential growth factors. 
This will cancel all but the exponential growth in each term. Since the upper bound and 
lower bound both have exponential growth factor (/? + 1), we get 

logos + 1) < lim - log^'(n) < logos + 1), 

n— >oo fl 



SO the exponential growth factor of q'{n) is OS + !)• 



n 
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So, we may now apply this to our example from before, getting a lower bound for the 
number of walks on Sy. 

Example 8.4. Take S5 and Sy, as we did in Example l8.3[ As discussed already, the number 
of walks on S7 obtained by inserting a north step at internal points of a walk from S5 is 

i>0 ^^' 
We showed in Section [8]X2l that jSj = 3, and it follows from qsin) < 3" that qsin) ~ C3"n'\ 
where C e R+ and s <0, therefore we may apply proposition |8.3.11| to find that 

q'j(n) ~ C'4V', 

where C must be positive real, since q'-j(n) counts something, and s' < since q'-j(n) < C4", 
as shown in Proposition l8.3.111 Thus, if qiin) is the number of all possible walks on S7, 
and h^{n) is the number of half plane walks taken on S7, we have 

q'jin) < q-,{n) < h-j{n). 

We may assume that n is sufficiently large, and replace each count with its asymptotic 
growth, giving 

CA"n'' < D/3'y < Ern". 

Using our standard method of isolating the exponential growth factor with a limit of a 
logarithm, we find 

log 4 < log;07 < log 4, 

or 

/St = 4. 

The previous example shows that the exponential growth factor for the sequence {qj{n))^ 
is 4. The next proposition gives us a method for adding two steps to arrive at a similar 
conclusion. 

Lemma 8.3.12. Let q{n) ~ Cp^n^ where C is a positive real constant, fi e [1, cc) and s <0. 
Let d(n) be the number of Dyck prefixes on n steps. Then the quantity 



q"(.n) = Y,hdii)q(n-i) 
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has growth 

q"in)-C"(J3 + 2)V", 

where C" is a positive real constant and s" < 0. 

Proof. Fix 6 > 0. Then thanks to the asymptotic form of q(n), we know there is some M > 
such that for all n > M,we have 

q(n) 



(1 - e) < 



Cj3"n' 



< (1 + e). 



For d(n), which has growth d(n) ~ D2"n'' with r < 0, there is some M' > such that for 

n > M' we have 

d(n) 



(1 - 6) < 



D2"n' 



< (1 + e). 



Take A^ = max{M', M} and let n > IN. Then we define 



q"{n) = Y^{l'\d{i)q{n-i). 

!>0 ^ ' 



The argument here follows the same lines as Proposition I8.3.5[ splitting the sum into 
three intervals, and bounding. Since it is similar, we move somewhat faster here. Going 
through the motions leads to the following inequality. 



(1-26) 



Cn' Yj \\d{i)P'-' - Dn'Xp'-') + CDn'^-'ifi + Tf 

+Dn' Yj \.Y^'q{n-i)-Cn'l'p'-') 

i=n-N+i ^ ' ) 



V i<N 



< q"(n). 



Now we analyse the growth. The argument from here follows the same lines as the 
single step case. We divide throughout by CD(J3 + ly^n'"^'^ to find 



(1-26) 



+ 1 



+ 



i:i.Nl;!)id{i)/3"-'-Dn'-2T-') 
D(J3 + lyn'- 

i:L-7..i0)(2'^(«-O-Cn^2';8"-0) 



C(J3 + lYn' 



< 



q"(n) 



CD(J3 + 2Yn''+' 
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The error terms both tend to zero, so we can find some A^~ large enough such that for n > N~ 
both will be smaller than e/2. So, taking n > N~ means 

CD(J3 + lYn'-^' 
The upper bound is similar, so again we omit some details. Since (1 + 3e) is larger 
in modulus then both (1 + 6) and (1 + 6)^, we replace both of the smaller terms by the 
larger, increasing the upper bound. We also replace the monomials in i and n by 1, since the 
exponents are at most 0. We then do the same analysis, allowing us to find some A^^ such 
that for n > A^+, 

(l+46)> 



q"{n) 



CD(j3 + 2)" 
Putting this all together, we find that for n > maxjA'^^, A'^'^} we have 

(1 - 4e)CD(J3 + Ifn'^' < q"{n) < (1 + Ae)CD{fi + If. 

Again, by using our method of isolating the exponential growth factor, we find that the 
exponential growth factor of q"{n) is jS + 2 as desired. The constant portion of the growth 
must be real and positive, since q"{n) is positive, and finally, the polynomial part must have 
exponent at most zero, since the upper bound has no polynomial portion. That is 

q"{n)'-C"{fi + 2fn'", 

where C" is positive and real and s" < 0. n 

Let us see this applied to an example. 

Example 8.5. Take again Sy and Sg as we did in Example 18.31 Now we wish to build 
walks of length n taken on the step set Sg by inserting a meander (or Dyck prefix) of length 
n - / on the {E, W} into a walk of length / on Sy. If q-j^i) is the number of walks of length / 
on Sy and d{n - i) the number of Dyck prefixes of length n - i, the number of such walks is 

Now that we have the proofs of the base cases and Lemmas 18.3.111 and 18.3.121 we know 
that the growth of qiin) is suitable for the application of Proposition 18.3. 121 which tells us 
that 

q^{n) '^ C 6 n , 
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where C" is real and positive, and s" is at most zero. 

Next, we know that q'^(n) forms a lower bound for q^in) 

q'lin) < q%{n) < h(n), 

since we can not get all walks of length n on Sg in this way. We take n sufficiently large, so 
that we may make the substitution of each quantity with its growth, giving 

C"6"n'" < DfiV < E&'n'. 

Once again, we may isolate the exponential growth factor, finding 

log 6 < logjS < log 6, 

or simply 

yS = 6. 

Thus, we have shown the exponential growth factor for the sequence (q^in)),, is fig, = 6. 

Thus we have a way of adding one step and two steps at a time to step sets and increasing 
the exponential growth. This allows us to give lower bounds on the growth of a lot of step 
sets, using fewer than half of the 23 D-finite cases as groundwork. Lemmas 18.3.111 and 
18.3. 12[ along with the results of Section I8.3.2[ give the results in Table 18. 1[ Before we 
conclude this chapter, we mention some other uses for our new lemmas. 

There is another way to use Lemma 18.3.121 that has not yet been discussed. So far, 
all examples show the insertion of a cardinality two step set with zero drift (steps in the 
opposite direction), but we can also insert two element step sets with positive drift, as long 
as no component of the drift is towards a boundary. Take S7 = i as an example. Rather 
than using the step set Y as a parent for this step set, we could use Lemma 18.3.1 2 I to insert a 
Dyck prefix on the steps \y into a Dyck prefix on the vertical steps I , as counting the ways 
to do this would give formula from the hypothesis of Lemma I8.3.12[ This gives a large 
subset of the walks on S7 as a shuffiing of Dyck prefixes (walks on Si), which combined 
with Theorem 18. 3. 101 gives the exponential growth factor of 4. 

This concludes the material of Chapter [8] and Part |ll] Part [III] will conclude the thesis 
with a discussion of the results and future work. 



Part III 



Conclusion 
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Chapter 9 

Summary and open problems 



The main result of this thesis is the proof of the exponential growth factors for the quarter 
plane models with either zero drift of a drift vector parallel to a boundary through combi- 
natorial means. This was done by bounding the number of walks on a step set S, by known 
quantities with exponential growth factors equal to the conjectured exponential growth of 
each model. By then taking a limit which isolated the exponential growth of each bound, 
we were able to squeeze the exponential growth and prove the conjectured value in each 
case. 

The significance of these results can be expressed in two ways. First, it answers any 
question about the possibility of a human proof of the exponential growth of quarter plane 
walks with single component drift in the positive. Secondly, up until now there has been 
no treatment of the exponential growth of walks taken on these 23 step sets in one piece 
of work, without either a level of computation similar to that given in ||4l[5l, or some very 
complicated analysis as in the forthcoming results announced by FayoUe and Raschel in 
[[Toll . Our method not only gives human driven, rigorous results using some fairly elemen- 
tary techniques, but it also introduces Lemmas 18.3.1 H and 18.3. 121 which give a systematic 
method for boot strapping exponential growth from smaller step sets. 

Our results are by no means comprehensive and for the remainder of this chapter we 
discuss extensions of these methods to more general settings. The first extension, and per- 
haps the lowest hanging fruit, is to work on the remaining 56 step sets - those with non-zero 
drift which is not parallel to an axis - using the method of bounding and squeezing devel- 
oped in Chapter[8] Next, we consider bounding subexponential growth factors, though also 
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Class 


Functional Equation 


Classification of GF 


Enumerative Results 


nVs 


Ws{x,y;t) = I + tS(x,y)Wzix,y;t) 


Rational 


Exact and asymptotic 
formulae for coefficients. 


Ts 


F^{y;t) = l + tP(y)F^(y;t) 
-{y'}tP(y)Fs{0; t) 


Algebraic 


Coefficient formulae from coefficient 
extraction of explicit GF, 


asymptotic results from Theorem 16.3.21 


■Ke 


Hdx,y;t) = \+tS(x,y)Hsix,y-t) 
-{y<'}tS(x,y)Hdx,0;t) 


Algebraic 


Results via bijection with !Fs 


Qe 


Qq(x, y;t)= l+tS (x, y)Qe(x, y; t) 

-t({x<'']S{x,y)Qr^{0,y;t) 

+{y^'']S{x,y)QQix,0-t)) 

+a;(S)^es(0,0;0 


23 D-finite 

5 1 non D-finite 

5 conjectured non D-finite 


For D-finite cases: 
enumerative results for 19 
cases via orbit sums, some direct 
counting results, experimental 
asymptotics via long sum 
generation, exponential growth 
factors from bounding arguments 



% 

Co 



g 
i 



Table 9.1: A table summarising each class, its functional equation, classification of GFs and summary of enumerative 
results. 



00 



CHAPTER 9. SUMMARY AND OPEN PROBLEMS 90 

address the prospect of limited success. Finally, we consider generalisations of Theorem 
18.3. 101 and Lemmas 18.3.1 II and l8.3. 121 to more general models: those with larger steps and 
those on higher dimensional lattices. 

Results for the remaining sets of small steps 

For the remaining sets of small steps, finding the exponential growth factors is harder. These 
models are known in many cases to have a non-D-finite generating function, meaning that 
they are not G-series. Therefore we don't have a nice template for the asymptotic growth of 
the counting sequence as in the 23 cases in Chapter [8] We can still apply the methodology 
as long as we assume a generic growth template 

qs{n) ~ K/3"e(n), 

where 6{n) is some subexponential factor, that is, lim„^oo \6{n)\^^" = 1. Again, here drift is 
the key in determining whether a model will have full exponential growth equal to the num- 
ber of steps, or exponential growth diminished by increased interaction with a boundary. 

Positive drift - tight bounds from Lemma 18.3.111 

We begin by showing a family of step sets known as the singular sets, shown in Figure 
19.11 for which our method will give exponential growth factors. Notice that all step sets 
have positive drift in each component, so the drift is into the quarter plane, and it would be 
reasonable to expect that in each case/?s = I Si- 





Figure 9.1: The singular step sets. 

These are all positive drift step sets, so the exponential growth for any counting se- 
quence is bounded above by the number of steps, found by taking the half plane walks 
bounded below by the jc-axis. We can bound the two step sets with cardinality 3 below by 
a subset of all walks, in both cases by injecting walks taken on step set Sn (or its y = x 
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reflection). Then since we know the results for the smallest two step sets, we can get the 
other three by applying Lemma [8.3.11l 

Negative drift - finding bounds requires care 

Now, we show an example S with more interesting drift, shown in Figure 19.21 The drift 
is positive in the y-component but negative in x. We can expect a diminished exponential 
growth factor due to this drift, but we must take care in bounding it. 




Figure 9.2: A non-singular step set S with a negative drift component. 

The drift is key to finding a least upper bound on the exponential growth. In Theorem 
18.3. 10[ it was shown that an upper bound on the exponential growth could be found by 
removing the y-axis as a boundary and using the results of Chapter |7]to find the exponential 
growth of the relaxation. What made these upper bounds the lowest available from this type 
of relaxation was that the drift was always parallel to the j-axis, and thus a negative drift 
quarter plane walk was also a negative drift half plane walk. 

However, as already noted, the drift is x-negative and y-positive. Thus, removing the 
^'-axis will give a positive drift half plane model, which won't give a least upper bound. 
Taking instead the relaxation to the x > half plane model, we preserve the negative drift 
and find the least upper bound available from this relaxation. 

Following the methods of Chapter[8l the next step would be to attempt to remove steps 
from S in order to find a more easily enumerated smaller model to which we could apply 
Lemma [8. 3. Ill or [8. 3. 121 There is only one choice of step to remove here, though, and it 
leads to the trivial model X . Injection of walks on Sn would once again work, but we 
cannot expect this to be tight. 

Some future work on this topic could focus on other ways of producing lower bounds 
to overcome the problems just discussed. One way was already touched on in the proof 
of Proposition I8.3.3[ where a lower bound on the jSj was imported from /3n by rotating a 
step away from a boundary. While this didn't help in the example just given, we could 
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find families related by this rotational relation to help cut down any case analysis needed. 
A long term goal would be finding an automatic way to bound growth constants below, 
as we did with upper bounds in Theorem 18.3. 101 Full asymptotic expressions may be on 
the way, however. Very recently, FayoUe and Raschel have announced 'general asymptotic 
results' in a forthcoming paper in [101 . giving as an example the asymptotic expression for 
the counting sequences of walks on Si ending anywhere, on an axis and returning to the 
origin. They do so by using integral equations found by viewing the model as a boundary 
value problem, requiring some heavy analysis to find solutions. 

Predicting polynomial growth 

Next, we move on to the prediction of the subexponential factors. To simplify the dis- 
cussion, we focus on our 23 models with simple drift. We know that for these cases, the 
subexponential factors are of the form 

e{n) = Krf, 

where /c is a positive real constant and a is, according to Bostan and Kauers, a negative 
rational number. For some models, such as S2 = {NE, SE, SW, NW}, we can find the 
polynomial factors by an application of Stirling's approximation, since we can count the 
whole set. For others, where we need to appeal to the asymptotics of a subset, we need a 
way to isolate the polynomial factor. The following is a very informal sequence of calcula- 
tions intended to show how we might do this. 

So, taking a step set S, we have our template for the asymptotics of the quarter plane 
walks on S 

q{n) ~ Kli"ff. 

For sufficiently large n we can replace q{n) with its asymptotic estimate and divide out the 
known exponential growth jS": 

— m \ 

/3" 

then taking a logarithm: 

logl — I -logK + alogn; 
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and the correct limit: 

lim ■; log — -— ~ a; 

n^oo log n \ /3" I 

we can isolate a. 

Given a method for isolation involving limits, it would be reasonable to expect to be 
able to bound and squeeze yet again. Unfortunately, the argument breaks down. Intuitively, 
it's precisely because the relaxations and restrictions have the same exponential growth: 
it's then the role of the polynomial factor to modulate the asymptotic growth and keep it 
close to the actual count. For the walks in the half plane, according to Theorem 16.3.21 if 
the drift is positive the polynomial factor is of degree 0. Similarly, the negative drift cases 
have polynomial factors of degree -|, and the zero drift have degree -^. In all cases, this is 
out by a difference of at least ^ with the results of Table 18.21 and will not give a tight upper 
bound on the polynomial factor. It appears, therefore, that new methods are necessary for 
finding tight bounds on the polynomial growth. Again, according to the announcement of 
FayoUe and Raschel in IfTOll . their approach using boundary value problems and integral 
equations will produce full asymptotic expressions for classes of walks with small steps. 

Generalisations 

The material of Chapter |6] is taken from a paper Ij2J, in which Banderier and Flajolet de- 
termine explicit enumerative formulae for directed paths with any step size. These more 
general results should give upper bounds on the exponential growth of the coefficients of 
generating functions for more general step sets in the quarter plane, meaning that Theorem 
18.3. 101 should translate directly. Furthermore, Lemmas 18.3.1 ll and l8. 3. 121 make no mention 
of the size of steps, and so they might be translated into a more general setting as well. This 
discussion leads to the following conjecture. 

Conjecture 1. Let S be as above. Let X = Q U {sq} and II = S U {^i, 52}, where sq is a 
step of any size away from any boundary and {^i, ^2} is a step set without drift towards any 
boundary. If ySs (respectively /3z,/3u) is the exponential growth of walks in the quarter plane 
on S (resp. Z, U), then Lemma l8.3. Ill will import from S the lower bound 
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and similarly. Lemma [8 . 3 . 1 2 1 will import the lower bound 

/33 + 2</3u. 

This conjecture is subject to the same considerations as its counterpart in Chapter|8] the 
application of Lemmas 18.3. Il l and 18.3. 121 correspond to the injection of single steps away 
from a boundary, or meanders on a pair of steps with drift away from a boundary. 

Another generalisation is the application of these results to higher dimensional lattice 
paths. This is an active area of research [^, with much experimental work being done 
at the time of writing. There should be no trouble in boot strapping as we have done for 
planar walks with small steps: as long as we are inserting either a single step away from any 
boundary or a pair of steps with either zero drift of non-zero drift away from any boundary, 
we should get analogous results. We summarise in the following conjecture. 

Conjecture 2. Let S be a set of small steps in Z'' for some dimension d, let R be some 
open region in Z'' bounded by rational hyperplanes, let J = S U {^o} and U = QU {si,S2} 
where sq is a step away from any boundary of R and si, sj is a pair of steps without drift 
towards any boundary. If ySe (respectively /3x,/3u) is the exponential growth of walks in R 
on S (resp. Z, U), then Lemma [8. 3. Ill will import from S the lower bound 

and similarly. Lemma 18.3.1 21 will import from S the lower bound 

/3e + 2</3ii. 

Finally, we have a conjecture about producing upper bounds. In the remaining quarter 
plane models, experimental work has shown that sometimes we must find other half planes - 
which contain the quarter plane - to produce upper bounds which are reasonable candidates 
for tight upper bounds. The conjecture is for step sets of any size in the quarter plane. 

Conjecture 3. If a step set S is symmetric about an axis L, with L passing through the first 
quadrant (potentially on a boundary), then the half plane walks bounded by the line L^ have 
the same exponential growth as the quarter plane walks. 
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